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Momentum  balances  on  the  North  Carolina  Inner  Shelf 

Steve  Lentz1,  R.  T.  Guza2,  S.  Elgar1,  F.  Feddersen2,  T.  H.  C.  Herbers3 

Abstract.  Four  months  of  moored  current,  pressure,  temperature,  conductivity, 
wave,  and  wind  observations  on  the  North  Carolina  shelf  indicate  three  dynamically 
distinct  regions:  the  surf  zone,  the  inner  shelf  between  the  surf  zone  and  the  13-m 
isobath,  and  the  mid  shelf.  In  the  surf  zone  the  along-shelf  momentum  balance 
is  between  the  cross-shelf  gradient  of  the  wave-radiation  stress  and  the  bottom 
stress.  The  linear  drag  coefficient  in  the  surf  zone  is  about  10  times  larger  than 
seawards  of  the  surf  zone.  On  the  inner  shelf  the  along-shelf  momentum  balance 
is  also  frictional;  the  along-shelf  wind  stress  and  pressure  gradient  are  balanced  by 
bottom  stress.  In  the  cross-shelf  momentum  balance,  the  pressure  gradient  is  the 
superposition  of  roughly  equal  contributions  from  the  Coriolis  force  (geostrophy) 
and  wave  setdown  from  shoaling,  unbroken  surface  gravity  waves.  At  mid  shelf 
the  along-shelf  momentum  balance  is  less  frictional  and  hence  flow  accelerations 
are  important.  The  cross-shelf  momentum  balance  is  predominantly  geostrophic 
because  the  greater  depth  and  smaller  bottom  slope  at  mid  shelf  reduces  the 
importance  of  wave  setdown.  The  cross-shelf  density  gradient  is  in  thermal  wind 
balance  with  the  vertical  shear  in  the  along-shelf  flow  in  depths  as  shallow  as 
10  m.  The  dominant  along-shelf  momentum  balances  provide  a  simple  estimate 
of  the  depth- averaged  along-shelf  current  in  terms  of  the  measured  forcing  (i.e. 
wind  stress,  wave  radiation  stress  divergence,  and  along-shelf  pressure  gradient) 
that  reproduces  accurately  the  observed  cross-shelf  variation  of  the  depth- averaged 
along-shelf  current  between  the  surf  zone  and  midshelf. 


1.  Introduction 

The  relative  importance  of  surface  gravity  wave  and  wind 
forcing  varies  by  an  order  of  magnitude  between  the  surf 
zone  (water  depths  of  order  1  m)  and  the  mid  shelf  (water 
depths  of  order  100  m).  Many  observational  programs  have 
focused  on  flows  in  the  surf  zone  or  on  the  mid  shelf.  How¬ 
ever,  there  are  few  detailed  observational  studies  of  the  inner 
shelf  region,  between  the  surf  zone  and  the  mid  shelf,  and 
momentum  balances  there  are  understood  poorly. 

The  depth-averaged  along-shelf  momentum  balance  within 
the  surf  zone  (assuming  along-shelf  homogeneous  bathymetry) 
is  between  forcing  by  obliquely  incident  breaking  surface 
waves  (e.g.  gradients  in  the  along-shelf  component  of  the 
wave  radiation  stress  dSxy/dx),  bottom  stress,  and  possi¬ 
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bly  cross-shelf  mixing  processes  [ Thornton  and  Guza,  1986; 
Svendsen  and Putrevu,  1994;  Feddersen  et  al.,  1998],  (Wave 
radiation  stresses,  e.g.  Sxy  and  Sxx,  represent  the  momen¬ 
tum  flux  due  to  surface  gravity  waves  and  are  analogous 
to  Reynold’s  stresses  [Longuet-Higgins  and  Stewart,  1964]. 
Unless  noted  otherwise,  the  time  scales  of  interest  are  a  few 
days,  long  compared  with  surface  wave  periods  and  vari¬ 
ables  are  averages  over  many  wave  periods).  Recent  studies 
[ Whitford  and  Thornton,  1993;  Feddersen  et  al.,  1998]  sug¬ 
gest  the  along-shelf  wind  stress  is  sometimes  significant  in 
the  surf  zone,  but  is  usually  much  smaller  than  dSxy/dx. 
The  along-shelf  momentum  balance  at  mid  shelf  is  complex. 
Acceleration,  Coriolis  force,  pressure  gradients,  wind  stress, 
and  bottom  stress  are  each  important  on  various  shelves 
[Allen  and  Smith,  1981;  Lentz  and  Winant,  1986;  Lee  et  al., 
1984,  1989].  Previous  studies  indicate  the  inner  shelf  is  dy¬ 
namically  different  with  along-shelf  flows  driven  by  both 
along-shelf  wind  stress  and  along-shelf  pressure  gradients, 
and  these  forcing  terms  are  balanced  primarily  by  bottom 
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friction  [Scott  and  Csanady,  1976;  Pettigrew ,  1981;  Lentz 
and  Winant,  1986;  Masse,  1988;  Lee  et  al.,  1989;  Lentz, 
1994], 

The  depth-averaged  cross-shelf  momentum  balance  within 
the  surf  zone  is  between  gradients  of  the  cross-shelf  radi¬ 
ation  stress  ( dSxx/dx )  associated  with  wave  breaking  and 
the  cross-shelf  pressure  gradient  dP/dx  ( i.e .  wave  set-up) 
[Longuet-Higgins  and  Stewart ,  1964;  Bowen  et  al.,  1968; 
Lentz  and  Raubenheimer ,  1999],  The  cross-shelf  momen¬ 
tum  balance  over  the  mid  shelf  is  between  dP/dx  and 
the  Coriolis  force  associated  with  the  along-shelf  flow  (i.e. 
geostrophic)  [Brown  etal.,  1985,  1987;  Lee  et  al.,  1989].  As 
the  depth  decreases  setup  (or  setdown)  forced  by  the  cross¬ 
shelf  wind  stress  will  become  increasingly  important.  On 
the  South  Carolina  shelf  in  10  m  of  water  dP/dx,  the  Cori¬ 
olis  force,  and  the  cross-shelf  wind  stress  were  estimated 
to  be  approximately  equal  in  magnitude  [Lee  et  al.,  1989]. 
However,  dP/dx  was  estimated  as  the  difference  between 
two  pressure  sensors  spanning  the  entire  75  km-wide  shelf. 
Furthermore,  surface  waves  were  not  measured  and  thus  the 
wave  setdown  associated  with  dSxx  / dx  in  unbroken  shoal¬ 
ing  waves  [Bowen  et  al.,  1968]  could  not  be  estimated.  The 
relative  importance  of  surface  wave  forcing,  Coriolis  force, 
and  cross-shelf  wind  stress  in  the  cross-shelf  momentum  bal¬ 
ance  over  the  inner  shelf  is  unknown. 

Previous  studies  have  generally  included  at  most  one 
mooring  site  between  the  5-m  and  30-m  isobaths  so  that  lit¬ 
tle  is  known  about  the  momentum  balances  across  the  in¬ 
ner  shelf.  More  detailed  observations  across  the  inner  16 
km  of  the  North  Carolina  shelf  are  used  here  to  determine 
the  cross-shelf  variation  of  the  dominant  terms  in  the  depth- 
averaged  cross-  and  along-shelf  momentum  balances.  The 
observations,  in  water  depths  ranging  from  4  to  26  m,  span 
the  region  from  the  surf  zone  to  mid  shelf  and  include  simul¬ 
taneous  measurements  of  currents,  winds,  waves,  and  pres¬ 
sures.  An  overview  of  the  observations  is  given  in  Section  2. 
Momentum  balances  are  described  in  Section  3.  Simplified 
scalings  of  the  observed  dominant  balances  are  considered 
in  Section  4,  followed  by  a  summary  in  Section  5. 

2.  Background 

2.1.  Field  Program  and  Data  Processing 

Observations  were  obtained  offshore  of  the  Army  Corps 
of  Engineers’  Field  Research  Facility  (FRF)  on  the  Outer 
Banks  near  Duck,  North  Carolina,  from  August  through 
early  December  1994  as  part  of  the  interdisciplinary  Coastal 
Ocean  Processes  Inner  Shelf  Study  [Butman,  1994].  The 
site  is  about  midway  between  Cape  Henry,  at  the  mouth  of 
Chesapeake  Bay  (100  km  to  the  north)  and  Cape  Hatteras. 
The  coastline  is  relatively  straight  between  Cape  Henry  and 
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Figure  1.  Plan  view  of  the  study  region  and  instrument  lo¬ 
cations.  The  central  cross-shelf  transect  of  instruments  is 
shown  in  detail  in  Figure  2. 

Oregon  Inlet  (50  km  to  the  south),  with  an  orientation  of 
about  340°T  at  the  array  center  (Figure  1).  The  bathymetry 
inshore  of  the  20-m  isobath  is  approximately  homogeneous 
along-shelf  on  scales  of  many  km.  Offshore  of  the  20-m  iso¬ 
bath  the  complex  ridge  and  swale  bathymetry  varies  several 
meters  vertically  over  horizontal  scales  of  a  few  km  (Figure 
2).  The  sea  floor  slopes  relatively  steeply  (0.01)  from  the 
4-m  to  the  13-m  isobath  and  more  gently  (0.002)  from  the 
13-m  to  the  20-m  isobath.  Offshore  of  the  20-m  isobath  the 
sea  floor  slope  is  small  (0.0002)  over  cross-shelf  scales  large 
compared  with  the  ridge  and  swale  features. 

Two  rigid  towers  (in  4  and  8  m  water  depth)  and  three 
surface/subsurface  mooring  pairs  (in  13,  21,  and  26  m  wa¬ 
ter  depth)  were  deployed  along  a  16-krn  cross-shelf  transect 
(Figure  2).  The  towers  supported  Marsh-McBirney  electro¬ 
magnetic  current  meters  and  fast-response  thermistors  sam¬ 
pled  at  2  Hz,  with  vertical  spacings  of  about  1  m  or  less.  The 
moorings  supported  Vector-Measuring  Current  Meters  that 
measure  horizontal  currents  and  temperature,  and  SeaBird 
SeaCATs  that  measure  temperature  and  conductivity,  all 
with  sample  rates  of  4  min.  Vertical  separations  for  these  in¬ 
struments  were  1  to  5  m  (Figure  2).  Two  SeaCATs  were  also 
deployed  1  and  4  m  above  the  bottom  on  a  piling  of  the  FRF 
pier  in  about  8-m  depth.  Setra  pressure  sensors,  sampled  at 
1  or  2  Hz  to  measure  surface  gravity  waves,  were  deployed 
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Figure  2.  Cross-shelf  instrument  transect  along  the  central 
line.  Thick  curve  is  the  bathymetry  (depth  relative  to  mean 
sea  level). 

about  1  m  above  the  sea  floor  near  each  mooring  or  tower, 
and  also  in  33-m  depth  (about  30-krn  offshore.  Figure  2). 
The  FRF  maintains  a  wave-directional  array  of  15  pressure 
sensors  in  8  m  depth  extending  about  200  m  along-shelf  and 
60  m  cross-shelf  and  sampled  at  2  Hz  [Long,  1996].  Wind 
velocity  was  measured  3.5  m  above  the  sea  surface  at  the 
21-m  site  (sample  rate  7.5  min)  and  19.5  m  above  the  sea 
surface  on  the  FRF  pier. 

Along-shelf  pressure  gradients  were  estimated  from  an 
array  of  eight  SeaBIRD  SeaGauges  that  each  measure  pres¬ 
sure,  temperature,  and  conductivity.  The  SeaGauges  were 
deployed  about  1  m  above  the  bottom,  with  five  along  the 
5-m  isobath  and  three  along  the  21-m  isobath,  in  both  cases 
spanning  an  along-shelf  distance  of  about  60  km  (Figure  1). 
Additionally,  at  the  northern  and  southern  21-m  sites,  Sea- 
CATs  were  mounted  1  m  below  the  sea  surface  on  surface 
buoys.  Along-shelf  array  instruments  were  sampled  every  4 
minutes.  Bottom  pressure  was  averaged  over  the  4  minute 
sample  interval.  Cross-shelf  pressure  gradients  on  the  array 
center  line  were  estimated  using  the  SeaGauges  on  the  5-m 
and  21-m  isobaths,  and  the  1  or  2  Hz  sampling  Setra  pres¬ 
sure  gauges  on  the  1 1-m,  26-m  and  33-m  isobaths  (Appendix 
A.l). 

The  instrumentation  was  deployed  during  late  July  and 
early  August  and  recovered  either  at  the  end  of  October 
or  in  early  December,  1994.  The  8-m  tower  and  the  13- 
m  surface  mooring  both  failed  about  October  10,  during  a 
strong  Nor’easter.  The  21-m  surface  mooring  failed  during 


a  Nor’easter  on  September  4,  came  ashore,  was  refurbished, 
and  redeployed  on  October  4.  Both  the  21-m  and  26-m  sur¬ 
face  moorings  failed  during  Hurricane  Gordon  (November 
16)  and  came  ashore.  The  instrumentation  and  data  were 
recovered  after  these  failures.  Time  series  from  the  4-m 
and  8-m  sites  have  gaps  owing  to  maintenance,  intermittent 
problems  with  individual  instruments,  and  a  severe  lightning 
strike  in  early  August. 

All  time  series  (with  the  exception  of  surface  wave  data) 
were  block-averaged  to  hourly  values  centered  on  the  hour. 
The  focus  is  on  subtidal  dynamics  so  time  series  were  low- 
pass  filtered  (half-power  point  38  hours)  unless  noted  oth¬ 
erwise.  All  vector  time  series  were  rotated  to  a  coordinate 
system  based  on  the  coastline  orientation  (Figure  1),  with 
the  along-shelf  coordinate  y  positive  toward  340°T,  and  the 
cross-shelf  coordinate  x  positive  offshore.  Obvious  biases 
and  drifts  in  some  conductivity  time  series  were  identified 
and  corrected  by  comparisons  with  adjacent  moored  con¬ 
ductivity  cells  and  shipboard  CTD  data  obtained  near  the 
moorings  [Alessi  et  al.,  1996;  Waldorf  et  al.,  1995,  1996]. 
Data  were  discarded  from  near-bottom  conductivity  cells 
that  drifted  substantially  in  late  October,  presumably  be¬ 
cause  of  fouling  by  suspended  sediment.  Salinity  and  den¬ 
sity  were  estimated  using  the  temperature  and  corrected  con¬ 
ductivity  [Fofonoff  and  Millard,  1983].  Pressure  time  series 
from  sensors  mounted  on  anchors  (water  depths  greater  than 
8  m)  sometimes  show  positive  shifts  of  1  -  30  mb  during 
storms,  presumably  owing  to  scouring  and  settling  of  the  an¬ 
chors.  Anchor  shifts  greater  than  1  mb  were  identified  and 
removed  by  comparison  with  time  series  from  other  near¬ 
bottom  pressure  sensors  (rigidly  mounted  on  jetted  pipes  in 
shallower  water)  that  did  not  shift.  A  more  detailed  descrip¬ 
tion  of  data  return  and  initial  processing,  including  correc¬ 
tion  of  conductivity  time  series  and  removal  of  anchor  shifts 
from  pressure  time  series,  is  given  in  Alessi  et  al.  [1996], 

2.2.  Overview  of  Observations 

The  dominant  time  scale  of  wind  variability  is  a  few  days, 
associated  primarily  with  the  passage  of  cold  fronts  [Austin, 
1998].  Comparison  of  winds  from  the  FRF  pier  (8-m  depth), 
the  21-m  site,  and  several  NDBC  buoys  in  the  region  indi¬ 
cate  that  subtidal  winds  did  not  vary  substantially  across  the 
16  km  cross-shelf  extent  of  the  study  region  during  the  field 
program  (Appendix  A. 2).  Wind  directions  are  typically  45° 
to  the  coast,  either  poleward  (up welling  favorable)  and  off¬ 
shore,  or  equatorward  (downwelling  favorable)  and  onshore. 
The  strongest  wind  stresses  (Appendix  A.l)  were  associated 
with  Nor’easters  (downwelling  favorable)  (Figure  3).  Sig¬ 
nificant  wave  heights  HSig  at  the  8-m  wave  array  ranged 
from  0.2  to  4  m  (Figure  3)  and  are  correlated  with  the  wind 
stress  magnitude.  Hourly  averaged  wave  directions  at  the  ar- 
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ray  ranged  between  +50°  and  —50°  from  normally  incident, 
i.e.  wave  crests  parallel  to  the  coast. 

Along-shelf  current  variance  over  time  scales  from  hours 
to  weeks  is  dominated  by  subtidal  variability.  At  all  five 
sites,  the  correlation  between  subtidal  along-shelf  currents  at 
different  levels  in  the  water  column  is  greater  than  0.75  (the 
95%  confidence  level  for  a  correlation  significantly  different 
from  zero  is  0.45  assuming  a  decorrelation  time  scale  of  3 
days  and  60-day  long  time  series).  Subtidal  along-shelf  cur¬ 
rents  tend  to  decrease  by  roughly  a  factor  of  two  from  near 
the  surface  to  near  the  bottom.  The  much  weaker  subtidal 
cross-shelf  currents  u  are  not  well  correlated  over  the  wa¬ 
ter  column  at  any  site,  and  are  often  in  opposite  directions 
near  the  surface  and  bottom.  Depth-averaged  and  subtidal 
(e.g.  low-pass  filtered.  Appendix  A.l)  currents  are  consid¬ 
ered  here  unless  otherwise  noted. 

The  mean,  depth-averaged,  along-shelf  flow  v  is  equator- 
ward  at  all  five  sites;  9  cm/s  at  the  13-m  site  and  4-5  cm/s  at 
the  4-m  and  26-m  sites  (Table  1).  Although  the  data  do  not 
span  the  same  time  period  (Figure  3),  the  cross-shelf  struc¬ 
ture  for  a  shorter  common  time  period  is  similar. 

The  principal  axes  of  the  depth-averaged  currents  are 
along-shelf  to  the  accuracy  of  the  measurements  («  5°) 
with  standard  deviations  of  13-20  cm/s  (Table  1).  Correla¬ 
tions  between  depth-averaged  along-shelf  currents  measured 
at  the  five  cross-shelf  sites  range  from  0.64  to  0.94.  In  con¬ 
trast  to  the  depth-averaged  along-shelf  flow,  both  the  means 
and  standard  deviations  of  the  depth-averaged  cross-shelf 
flow  u  are  small  relative  to  the  accuracy  of  the  current  mea¬ 
surements  (2-3  cm/s)  [Beardsley,  1987;  Guza  el  al.,  1988]. 
Depth-averaged  cross-shelf  currents  at  different  moorings 
are  generally  not  correlated,  consistent  with  error-dominated 
measurements  and/or  unresolved  spatial  variation  of  the  cur¬ 
rent.  In  either  case,  estimates  of  terms  in  the  momentum 
balances  associated  with  u  must  be  interpreted  cautiously. 
Though  a  few  cm/s  or  less,  u  may  still  be  important  for  cross¬ 
shelf  exchange  because  the  inner  shelf  region  is  so  narrow. 

The  vertical  structure  of  temperature  evolved  dramati¬ 
cally  during  the  study  (Figure  4a).  In  August,  23°C  near¬ 
surface  water  and  deeper  17°C  water  were  separated  by  a 
strong  thermocline  centered  about  10  m  below  the  surface. 
Wind-driven  upwelling  (downwelling)  of  the  thermocline  in 
August  resulted  in  large  cross-shelf  temperature  gradients 
onshore  of  the  26-m  site  as  the  thermocline  shoaled  (deep¬ 
ened)  and  formed  a  surface  (bottom)  front.  The  water  col¬ 
umn  was  vertically  mixed  at  least  as  far  offshore  as  the  26- 
m  site  in  response  to  strong  winds  from  the  northeast  in 
early  September.  A  strong  thermocline  did  not  redevelop 
and  temperature  differences  across  the  water  column  often 
were  much  less  than  2°C  from  September  through  Novem¬ 
ber. 


Despite  the  strong  thermocline  in  August,  both  vertical 
(Figure  4)  and  cross-shelf  density  gradients  were  dominated 
by  salinity  variations.  The  primary  source  of  salinity  vari¬ 
ability  was  narrow,  shallow  plumes  of  relatively  fresh  water, 
presumably  from  Chesapeake  Bay  [Rennie  el  al.,  1999].  The 
Chesapeake  Bay  plume  flowed  into  the  study  region  gener¬ 
ally  during  downwelling  favorable  winds  and  was  typically 
confined  inshore  of  the  21-m  isobath  until  it  was  swept  off¬ 
shore  by  upwelling  favorable  winds.  When  the  Chesapeake 
Bay  plume  was  not  present,  salinities  generally  increase  with 
depth  and  distance  offshore,  consistent  with  historical  data 
[Boicourt,  1973].  Vertical  salinity  (and  hence  density)  gra¬ 
dients  tended  to  be  large  during  August  and  small  during 
October  when  there  was  strong  wind  and  wave  forcing  (Fig¬ 
ure  3). 


3.  Momentum  Balances 


The  depth-averaged  momentum  equations,  assuming  hy¬ 
drostatic  flow  and  small  sea  level  variations  compared  with 
the  water  depth,  are 


du 

dt 


id  r°  2j  id  r°  , 

+  he-J.kndz  +  ml.rdz~lv 

_  1  dP  Tsx  _  rbx  _  1  dSxx 
p0  dx  p0h  pah  pah  dx 


dv 

dt 


+ 


1  d 
h  dx 


uvdz  + 


Id 

hdy 


v2dz  +  fu 


j_dP_  t^_  _  t^_  _  i  dsxy 

Po  dy  poh  Poh  Poh  dx 


(2) 


where  (u,v)  are  the  cross-shelf  (x)  and  along-shelf  (y)  sub¬ 
tidal  components  of  velocity,  ( u,v )  are  the  correspond¬ 
ing  depth-averaged  velocities,  2  is  height  above  mean  sea 
level,  h  is  the  water  depth,  /  =  8.59  x  10-5  s-1  is  the 
Coriolis  parameter,  pa  =  1023  kg/m3  is  a  reference  den¬ 
sity,  ( dP/dx ,  dP/dy)  is  the  depth-averaged  horizontal  pres¬ 
sure  gradient,  g  =  9.81  m/s2  is  gravitational  acceleration, 
(■ tsx ,  rsy)  is  the  wind  stress,  (rbx  ,rby)  is  the  bottom  stress, 
and  Sxx  and  Sxy  are  wave  radiation  stresses.  Lateral  mixing 
processes  on  time  scales  shorter  than  subtidal  may  be  signifi¬ 
cant,  especially  in  the  surf  zone,  but  are  not  included  because 
they  could  not  be  estimated  accurately  from  the  observa¬ 
tions.  Along-shelf  gradients  in  the  surf  zone  radiation  stress 
and  pressure  fields  at  relatively  short  spatial  scales  (e.g.  as¬ 
sociated  with  local  along-shelf  bathymetric  variations)  could 
not  be  estimated.  Support  for  neglecting  along-shelf  varia¬ 
tions  in  the  along-shelf  momentum  balance  (2)  inshore  of 
the  8-m  isobath  is  provided  by  Feddersen  el  al.  [1998], 


along-shelf  current  cm/s  N/m 
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Figure  3.  Time  series  of  the  along-shelf  wind  stress  (negative  is  equatorward),  significant  wave  height  in  8  m  water  depth, 
and  the  subtidal  depth-averaged  along-shelf  current  at  the  five  mooring  sites.  Offshore  distances  of  mooring  sites  are  in 
parentheses. 
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Site 

Mean 

Principal  Axes 

Days 

u 

V 

Major 

Minor 

Orientation,  deg 

4  m 

2 

-4 

18 

6 

4° 

79 

8  m 

1 

-7 

20 

2 

-1° 

61 

13  m 

-1 

-9 

19 

2 

-4° 

66 

21  m 

-1 

-8 

18 

2 

-3° 

69 

26  m 

-2 

-5 

13 

2 

2° 

93 

Table  1.  Statistics  of  depth-averaged  currents  (cm/s).  Orientations  are  relative  to  the  along-shelf  direction  340°T.  Time  series 
span  different  time  periods  (see  Figure  3). 


a)  temperature  (°C) 


b)  salinity  (psu) 


with  a  corresponding  equation  for  dPj d y.  dPs /dx  includes 
contributions  from  the  sea  surface  slope  and  the  atmospheric 
pressure  gradient  and  will  be  referred  to  as  the  barotropic 
pressure  gradient.  The  second  term  on  the  right-hand  side  of 
(3),  the  density  contribution  to  the  depth-averaged  pressure 
gradient,  will  be  referred  to  as  the  baroclinic  pressure  gradi¬ 
ent.  Assuming  hydrostatic  flow,  the  near-surface  pressure  in 
(3)  may  be  expressed  in  terms  of  near-bottom  pressure  Pb 
and  density 

,° 

ps  _  pb  _  /  pgdz.  (4) 

J-h 

The  terms  in  (1-3)  are  estimated  at  some  or  all  of  the 
five  mooring  sites.  Cross-shelf  gradients  are  estimated  as 
finite  differences  centered  on  the  mooring  sites.  An  alter¬ 
nate  approach,  estimating  terms  in  control  volumes  bounded 
by  adjacent  moorings,  yields  similar  results.  The  nonlin¬ 
ear  terms  on  the  left-hand  side  of  (1)  and  (2)  could  not  be 
estimated  accurately  at  the  4-m  and  8-m  sites  because  the 
wave  height  can  be  a  significant  fraction  of  the  water  depth 
and  the  water  column  above  the  wave  trough  level  (where  u 
is  onshore)  was  not  sampled.  Crude  estimates  suggest  the 
nonlinear  terms  at  the  deeper  sites  were  small.  The  non¬ 
linear  terms  are  not  considered  further.  Estimation  of  the 
other  terms  is  described  in  Appendix  A.l  and  uncertainties 
in  some  estimates  are  discussed  in  Appendix  A. 2.  Bottom 
stress  is  estimated  using  a  linear  drag  law 


Figure  4.  Time  series  of  (a)  temperature,  (b)  salinity,  and  (c) 
oo  for  depths  below  the  surface  of  2  m,  14  m,  and  24.5  m  at 
the  26-m  site. 

The  depth-averaged  pressure  gradient  (obtained  by  verti¬ 
cally  integrating  the  hydrostatic  equation)  includes  surface 
pressure  Ps  and  density  p  contributions 


(Tbx,Tbv)  =  p0r(ub,vb)  (5) 

where  r  —  5  x  10-4  m/s  based  on  previous  mid-shelf  ob¬ 
servations  [e.g.  Lentz  and  Winant ,  1986]  and  (ub,vb)  is  the 
velocity  0.5  to  1.5  m  above  the  bottom.  Bottom  stress  esti¬ 
mates  from  (5)  are  similar  to  log-profile  estimates  using  bot¬ 
tom  tripod  measurements  at  the  21-m  site  (Appendix  A. 3), 
suggesting  a  constant  r  provides  reasonable  bottom  stress 
estimates  at  this  site.  However,  the  along-shelf  momentum 
balance  discussed  in  Section  3.1  implies  r  is  substantially 
larger  in  the  surf  zone. 


LENTZ  ETAL.:  MOMENTUM  BALANCES  ON  THE  NORTH  CAROLINA  INNER  SHELF 


7 


The  time-average  over  the  deployment  period  of  terms  in 
the  cross-  and  along-shelf  momentum  balances  either  could 
not  be  estimated  ( e.g .  mean  pressure  gradients)  or  are  small 
compared  with  the  fluctuations  in  those  terms  and  are  not 
considered  further. 

3.1.  Along-shelf  Momentum  Balance 

At  the  4-m  site,  time  series  of  the  along-shelf  wind  ( rsy ) 
and  bottom  stresses  ( rby  from  (5))  are  similar  (Figure  5a). 
Exceptions  occur  during  portions  of  October  14-15  and 
September  22  when  the  along-shelf  wind  and  along-shelf 
current  are  opposed  (e.g.  Tsy  and  rby  have  opposite  signs 
in  Figure  5a,  and  in  Figure  8a  discussed  below).  At  these 
times,  the  wave  radiation  stress  gradient  opposed  the  wind 
stress,  the  4-m  site  was  in  the  surf  zone,  and  the  4-m  along- 
shelf  current  was  in  the  direction  of  the  wave  radiation  stress 
gradient  indicating  the  dominant  driving  force  was  the  radi¬ 
ation  stress  gradient  dSxy/dx  resulting  from  obliquely  in¬ 
cident  breaking  waves  (Figure  5b).  To  examine  where  wind 
or  waves  dominate  the  along-shelf  forcing  it  is  convenient  to 
categorize  each  site  (for  each  hourly  time  period)  as  either 
seawards  of,  or  in,  the  surf  zone,  even  though  the  demarca¬ 
tion  between  the  surf  zone  and  the  inner  shelf  is  not  shaip. 
For  simplicity,  a  site  is  considered  seawards  of  the  surf  zone 
if  HSig/h  <  0.33,  corresponding  to  breaking  of  less  than 
about  10%  of  the  waves  [Thornton  and  Guza,  1984].  Surface 
wave  heights  at  each  site  were  estimated  using  a  nearby  bot¬ 
tom  pressure  gauge  and  linear  theory.  Using  this  criterion, 
the  4-m  site  was  in  the  surf  zone  on  August  23,  September 
3-6,  19,  22,  October  3-4,  and  10-20.  The  4-m  site  is  within 
the  surf  zone  when  dSxy/dx  is  large  (Figure  5b)  and  also 
when  iTgjg  is  large,  but  the  wave  direction  is  close  to  nor¬ 
mally  incident  so  dSxy/dx  (A3)  is  small.  The  8-m  tower 
was  in  the  surf  zone  for  a  few  hours  on  September  4  and  22. 
The  21-m  and  26-m  sites  were  never  in  the  surf  zone,  and 
the  8-m  tower  and  13-m  surface  mooring  failed  during  the 
large  waves  in  mid  October.  More  complicated  expressions 
for  defining  the  offshore  extent  of  the  surf  zone  [e.g.  Bat- 
tjes  and  Stive ,  1985]  yield  similar  estimates,  and  the  results 
below  do  not  depend  critically  on  this  definition. 

When  the  4-m  site  is  in  the  surf  zone  the  standard  de¬ 
viation  of  (poh)-1  dSxy/ dx  is  1-2  orders  of  magnitude 
larger  than  the  standard  deviations  of  other  terms  in  the 
along-shelf  momentum  balance  (Table  2,  Figure  5).  To  bal¬ 
ance  dSxy/dx  with  the  bottom  stress,  the  linear  drag  co¬ 
efficient  r  for  the  subtidal  flow  in  the  surf  zone  must  be 
r  =  5  x  10-3  m/s,  10  times  larger  than  the  nominal  value  of 
r  (Figure  6).  The  linear  drag  coefficient  could  increase  be¬ 
cause  breaking  waves  transfer  momentum  effectively  to  the 
bottom,  or  because  bottom  roughness  increases  inside  the 
surf  zone  [Garcez-Faria  et  al.,  1998].  The  increase  in  the 


Site 

dv 

dt 

fu 

1  dP 
Po  dy 

Tsy 

Po  h 

rby 

Po  h 

1  dSvy 

p0h  dx 

4-m 

0.2 

0.2 

0.3 

0.9 

1.2 

0.1  (12.9) 

8-m 

0.3 

0.1 

- 

0.8 

0.7 

0.1  (0.1) 

13-m 

0.3 

0.2 

- 

0.4 

0.4 

0.0 

21-m 

0.2 

0.2 

0.3 

0.4 

0.3 

0.0 

26-m 

0.2 

0.2 

- 

0.3 

0.2 

0.0 

Table  2.  Standard  deviations  of  terms  in  the  along-shelf  mo¬ 
mentum  balance  (2).  Statistics  for  the  4-m  and  8-m  sites  are 
for  periods  when  these  sites  are  seawards  of  the  surf  zone. 
The  values  in  parentheses  for  ( p0h)~ldSxy/dx  include  all 
periods.  Units  are  10-5  m/s2. 

drag  coefficient  in  the  surf  zone  is  qualitatively  consistent 
with  results  of  a  more  detailed  examination  of  the  along- 
shelf  momentum  balance  onshore  of  the  8-m  site  using  addi¬ 
tional  observations,  a  cross-shelf  integration  that  reduces  the 
effects  of  lateral  mixing  on  the  estimated  drag  coefficients, 
and  a  quadratic  drag  formulation  that  includes  the  effect  of 
surface  waves  on  the  mean  bottom  stress  [Feddersen  et  al., 
1998].  Seaward  of  the  surf  zone  the  estimated  dSxy/dx  term 
is  small  (Table  2). 

The  estimated  Coriolis  force  fu  is  relatively  small  at  the 
4-m  and  8-m  sites,  but  similar  to  the  other  terms  at  the  21- 
m  and  26-m  sites  (Table  2).  However,  fu  will  be  neglected 
in  subsequent  analysis  because  estimates  of  fu  are  smaller 
than  the  estimated  uncertainties  (Appendix  A.2),  and  un¬ 
correlated  between  sites.  Furthermore,  fu  is  not  correlated 
(range  0.09  to  0.37)  with  the  the  sum  of  the  acceleration, 
pressure  gradient,  and  surface  and  bottom  stresses.  In  con¬ 
trast,  the  acceleration  term,  which  has  similar  magnitudes, 
is  correlated  (range  0.55  to  0.80)  with  the  sum  of  the  Corio¬ 
lis  and  surface  and  bottom  stress  terms  at  all  but  the  4-m  site 
(correlation  0.36).  These  results  are  consistent  with  previous 
studies  [Allen  and  Smith,  1981;  Lentz  and  Winant,  1986;  Pet¬ 
tigrew,  1981;  Lentz,  1994]  and  suggest  fu  is  not  estimated 
accurately  from  the  observations. 

When  the  4-m  site  is  offshore  of  the  surf  zone  the  stan¬ 
dard  deviations  of  the  surface  wind  and  bottom  stresses  are 
at  least  three  times  as  large  as  the  other  terms  (Table  2)  and 
correlated  with  a  regression  coefficient  near  1.0  (Table  3). 
Time  series  of  surface  and  bottom  stress  are  correlated  at 
all  sites  (Table  3).  However,  standard  deviations  of  dv/dt 
and  dP/dy  become  increasingly  significant  as  the  depth  in¬ 
creases  and  tsv / ( p0h )  and  rby /  (p0h)  decrease  (Table  2).  At 
the  8-m  site  the  standard  deviation  of  the  along-shelf  flow 
acceleration  and  pressure  gradient  (measured  along  the  5-m 
isobath)  are  similar  in  magnitude,  but  are  only  about  half 
the  magnitude  of  the  surface  and  bottom  stress  terms  (Ta- 
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a)  Tsy/psh  and  iby/p„h  4-m  site 


Figure  5.  Time  series  (unfiltered  hourly  values)  of  largest  terms  in  the  along-shelf  momentum  balance  at  the  4-m  site,  (a)  wind 
tsv / p0h  and  bottom  rbv / p0h  stresses,  and  (b)  radiation  stress  gradient  —  (p0h)~1dSxy/dx.  When  \dSxy/dx\  is  relatively 
large  the  site  is  within  the  surf  zone.  Note  the  vertical  scales  in  (a)  and  (b)  differ. 


ble  2).  The  dP/dy  term  along  the  21-m  isobath,  and  the 
acceleration  and  stress  terms  at  the  13-m,  21-m,  and  26-m 
sites  are  all  of  similar  magnitude  (Table  2,  Figure  7).  When 
the  along-shelf  pressure  gradient  is  largest,  the  surface  and 
bottom  stresses  differ  substantially  (mid  August,  September 
5,  October  16  in  Figure  7).  The  salinity  observations  indicate 
the  positive  dP/dy  events  in  August  that  drive  a  southeast¬ 
ward  flow  are  associated  with  the  Chesapeake  Bay  plume 
[Rennie  el  al.,  1999],  The  along-shelf  momentum  balance 
approximately  closes  at  each  site.  The  sum  of  the  response 
terms  (dv/dt  +  (p0h)~1Tby)  is  well  correlated  with  the  sum 
of  the  two  forcing  terms  (r?p„'  1  dP/dy  +  ( p0h)~  tsv)  at 
each  site,  and  the  regression  coefficients  are  about  1 .0  (Fig¬ 
ure  8,  Table  3).  Furthermore,  at  each  site  these  correlations 
are  higher  than  the  correlations  between  the  surface  and  bot¬ 
tom  stress  alone. 


Along  the  5-m  isobath  the  standard  deviation  of  the  barotropic 
pressure  gradient  is  about  five  times  larger  than  the  standard 
deviation  of  the  baroclinic  component.  In  contrast,  along 
the  21-m  isobath  the  standard  deviation  of  the  barotropic 
pressure  gradient  is  about  twice  the  baroclinic  component 
and  the  baroclinic  component  tends  to  oppose  the  barotropic 
component.  The  baroclinic  pressure  gradient  results  primar¬ 
ily  from  salinity  variations  and  is  substantial  even  when  the 
water  column  is  not  stratified  ( e.g .  early  October).  These  es¬ 
timates  indicate  that,  at  least  along  the  21-m  isobath  at  this 
site,  density  estimates  are  critical  for  accurate  estimation  of 
along-shelf  pressure  gradients. 
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Tby/(p0h)  vs  rsy 

Up  oh) 

R  vs  F 

Site 

a 

b 

Correlation 

a 

b 

Correlation 

4-m 

0.94±0.46 

0.08±0.41 

0.69 

0.98±0.36 

0.01±0.38 

0.82 

8-m 

0.71  ±0.25 

-0.13±0.20 

0.78 

0.86±0.20 

-0.02±0.18 

0.89 

13-m 

0.63±0.18 

-0.06±0.10 

0.80 

0.88±0.17 

0.02±0.10 

0.91 

21-m 

0.54±0.16 

0.01  ±0.06 

0.72 

0.75±0.17 

0.08±0.08 

0.88 

26-m 

0.56±0.17 

0.01  ±0.05 

0.73 

0.78±0.15 

0.08±0.05 

0.87 

Table  3.  Results  of  linear  regression  analysis  of  forcing,  F  =  (—dP/dy)/p0  ±  rsy /(p0h),  with  response,  R  =  dv/dt,  ± 
rby/(Poh),  in  along-shelf  momentum  balance  of  the  form  R  =  uF  ±  b.  Results  of  regression  between  Tsy/(p0h)  and 
Tby/(p0h)  are  also  given.  Analysis  includes  only  periods  when  sites  are  seawards  of  the  surf  zone.  All  correlations  are 
different  from  zero  at  the  95%  confidence  level  and  95%  confidence  for  a  and  b  are  shown.  Units  are  10-5m/s2  for  intercepts. 


Site 

du 

dt 

fv 

1  dP 

po  dx 

sx 

T _ 

Po  h 

TbX 

Po  h 

1  dSxx 

p0h  dx 

4-m 

0.0 

1.1 

- 

0.9 

0.4 

9.6(221.5) 

8-m 

0.0 

1.7 

- 

0.8 

0.2 

3.6  (3.5) 

13-m 

0.0 

1.6 

2.1 

0.4 

0.1 

1.5 

21-m 

0.0 

1.4 

1.8 

0.3 

0.1 

0.2 

26-m 

0.1 

1.1 

1.4 

0.3 

0.1 

0.6 

Table  4.  Standard  deviations  of  terms  in  the  cross-shelf  momentum  balance  (1).  Statistics  for  the  4-m  and  8-m  sites  are  for 
periods  when  these  sites  are  seawards  of  the  surf  zone.  The  values  in  parentheses  for  (pah~1)dSxx/dx  include  all  periods. 
Units  are  10_5m/s2. 


3.2.  Cross-shelf  Momentum  Balance 

In  the  surf  zone  the  standard  deviation  of  (pa  h )  1  dSxx/ dx 
is  several  orders  of  magnitude  larger  than  standard  devia¬ 
tions  of  the  other  terms  in  the  cross-shelf  momentum  bal¬ 
ance  (Table  4).  Presumably  dSxx/dx  is  balanced  by  a  cross¬ 
shelf  pressure  gradient  dP/dx,  consistent  with  wave-driven 
setup.  This  balance  cannot  be  confirmed  because  pressure 
gradients  were  not  measured  accurately  in  this  region  (Ap¬ 
pendix  A. 2).  However,  analysis  of  a  3.5  year  long  time  se¬ 
ries  (including  the  period  of  this  study)  of  pressure  in  2  m 
and  8  m  depth  at  the  Duck  site  indicates  that  large  negative 
9Sxx/dx  across  the  surf  zone  are  balanced  by  dP/dx,  con¬ 
sistent  with  wave  setup  [Lentz  and Raubenheimer,  1999]. 

When  the  4-m  site  is  offshore  of  the  surf  zone,  the  stan¬ 
dard  deviation  of  ( p0h)~1dSxx/dx  is  an  order  of  magni¬ 
tude  smaller  than  when  the  site  is  within  the  surf  zone,  but 
still  an  order  of  magnitude  larger  than  the  standard  devia¬ 
tions  of  the  other  estimated  terms,  suggesting  that  dP/dx 
is  dominated  by  (unmeasured)  wave  setdown.  At  the  8- 
m  site  the  magnitude  of  (/j0/i)_1  dSxx  / dx  is  a  factor  of 
three  larger  than  fv  and  the  cross-shelf  wind  stress  term 
(p0h)~1Tsx  (Table  4  and  Figure  9).  The  Coriolis  force  fv 
and  the  cross-shelf  wind  stress  tsx  are  positively  correlated 


(Figure  9a).  The  strongest  winds  are  from  the  northeast  and 
drive  an  equatorward  along-shelf  current  with  an  onshore  di¬ 
rected  Coriolis  force  that  reinforces  the  onshore  component 
of  the  wind  stress  so  that  the  largest  magnitudes  of  fv  and 
(/90(i)_1rsx  are  negative  and  occur  concurrently  (Figure  9a). 
However,  the  8-m  site  is  usually  seawards  of  the  surf  zone 
and  the  wave-driven  setdown  (dSxx/dx  is  negative)  opposes 
and  is  partially  balanced  by  the  wind  and  Coriolis  forced 
setup  during  Nor’easters  (Figure  9b).  Pressure  gradients 
-Po^dP/dx  =  ~(p0h)~1dSxx/dx  +  fv  +  ( p0h)~lrsx 
during  these  events  should  be  smaller  than  those  from  wave 
setdown  alone,  and  should  have  sign  opposite  to  those  from 
fv  +  (p0/i)_1Tsx  alone.  Again,  accurate  estimates  of  the 
cross-shelf  pressure  gradient  are  needed  to  confirm  this  bal¬ 
ance.  Inside  the  surf  zone  dSxx/dx  is  positive,  and  during  a 
Nor’easter  wave  setup  is  reinforced  by  the  cross-shelf  wind 
stress  and  Coriolis  force.  However,  the  dSxx/dx  term  is  so 
much  larger  than  the  other  terms  that  they  are  likely  negligi¬ 
ble  for  most  surf  zone  conditions  (see  Section  4.2). 

The  cross-shelf  pressure  gradient  dP/dx  was  measured 
at  the  13-m  site,  and  is  balanced  by  ( p0h)~  dSxx/dx  and 
fv.  The  cross-shelf  wind  stress  term  tsx  (p0h)~  is  reduced 
by  about  40%  at  the  13-m  site  relative  to  the  8-m  site  due 
to  the  increase  in  water  depth  (Figure  10a,  Table  4).  As  in 
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4-m  site 


-10  -5  0  .5  10  15 

-(3Sxy/ax)/p0h  (1 0-4  m/s2) 


Figure  6.  Hourly  values  of  along-shelf  stresses  Tby  / p0  h  ver¬ 
sus  —  (poJTj-hdSxy/dx  when  the  4-m  site  is  within  the  surf 
zone.  The  linear  relationship  and  the  slope  of  0. 1  suggests 
the  linear  drag  coefficient  r  within  the  surf  zone  is  10  times 
larger  than  seawards  of  the  surf  zone. 


Site 

Correlation 

Coefficient 

Regression 

Coefficient 

Days 

13-m 

0.63 

0.82  ±  0.38 

63 

1 13-m 

0.80 

1.34  ±  0.41 

48 

21-m 

0.92 

1.16  ±  0.22 

50 

26-m 

0.93 

1.19  ±0.15 

78 

Table  5.  Regression  analysis  of  fv  with  (dP/dx)  / p„, 
terms  in  cross-shelf  geostrophic  balance.  All  correlations 
are  significantly  different  from  zero  at  the  95%  confidence 
level  and  95%  confidence  intervals  for  regression  coeffi¬ 
cients  are  shown,  f  Regression  analysis  of  fv  —  with 
(dP/dx)/ pa. 

8-m  depth,  —  (pohfj^dSxx/dx  and  fv  +  (p0/i)_1Tsx  are 
opposed  during  the  strongest  events,  and  in  13-m  depth  can¬ 
cellation  is  nearly  complete  in  some  cases  ( e.g .  September  5, 
October  4  and  October  10-11  in  Figure  10b).  The  observed 
pressure  gradients  are  qualitatively  consistent  with  this  can¬ 
cellation  (Figure  10c).  The  Coriolis  force  fv,  dSxx/dx, 
and  dP/dx  are  mutually  correlated.  Linear  regression  of 
fv  and  pa~1dP  /  dx  yields  a  significant  correlation  (0.63) 
and  a  regression  coefficient  of  0.82  (Table  5),  suggesting  a 
geostrophic  balance.  However,  pa~ldP  /  dx  is  better  cor¬ 
related  (0.80)  with  (fv  —  ( p0h)~ldSxx/dx )  than  with  fv 
alone.  These  results  (Tables  4  and  5,  and  Figures  9  and  10) 


suggest  that  8Sxx/dx  is  important  in  the  cross-shelf  mo¬ 
mentum  balance  offshore  of  the  surf  zone,  in  depths  at  least 
as  great  as  13  m. 

At  the  21-m  and  26-m  sites  fv  and  pa~ldP/dx  are 
highly  correlated  with  linear  regression  coefficients  near  one 
(Table  5)  and  are  larger  than  the  other  terms  (Table  4). 
The  cross-shelf  momentum  balance  is  thus  consistent  with 
geostrophy  (Figure  11).  The  standard  deviations  of  rsx / pah 
and  ( p0h)~  dSxx/dx  are  less  than  1/5  of  fv  in  21-m  depth 
(Table  4),  and  the  correlation  between  fv  and  p0~ldP  /  dx 
is  not  altered  significantly  by  including  either  in  the  bal¬ 
ance.  At  the  26-m  site,  the  standard  deviation  of  the  esti¬ 
mated  ( p0h)~1dSxx/dx  is  larger  (about  1/2  of  fv)  than  at 
the  21-m  site  because  the  local  bottom  slope  is  steeper  (hx 
in  (15),  discussed  below).  However,  the  correlation  is  not 
increased  by  including  dSxx/dx,  possibly  because  of  the 
disparity  in  the  spatial  scales  used  to  estimate  dSxx/9x  and 
dP/dx.  The  local  bottom  slope  in  26  m  depth  is  used  to 
estimate  dSxx/dx,  but  the  local  slope  is  larger  and  of  differ¬ 
ent  sign  than  the  average  bottom  slope  over  the  cross-shelf 
separation  of  the  pressure  gauges  (at  the  5-m  and  33-m  sites) 
used  to  estimate  dP/dx  at  the  26-m  site  (Figure  2). 

The  baroclinic  component  of  the  cross-shelf  pressure  gra¬ 
dient  often  opposes  the  barotropic  component  at  the  13-m 
and  21-m  sites  where  bottom  pressure  gradient  and  den¬ 
sity  gradient  estimates  could  be  made  (Figure  12).  Corre¬ 
lations  between  the  barotropic  and  baroclinic  pressure  gra¬ 
dients  are  -0.82  at  both  sites.  The  standard  deviations  in¬ 
dicate  that  the  depth-averaged  baroclinic  pressure  gradient 
balances  30-45%  of  the  barotropic  pressure  gradient.  Cross¬ 
shelf  density  gradients  (baroclinic  pressure  gradients)  in  Oc¬ 
tober,  when  stratification  is  relatively  weak  (Figure  4),  are 
similar  in  magnitude  to  cross-shelf  density  gradients  in  Au¬ 
gust  when  stratification  is  relatively  strong  (Figure  12). 

The  importance  of  the  baroclinic  contribution  to  the  depth- 
averaged  cross-shelf  pressure  gradient  and  the  tendency  for  a 
geostrophic  balance  suggests  the  vertical  shear  in  the  along- 
shelf  velocity  may  be  in  thermal  wind  balance  over  the  inner 
shelf, 

9v_  = _ 9_dp_ 

dz  p0f  dx' 

Although  (6)  has  been  shown  to  hold  in  the  middle  of  the 
water  column  at  mid  shelf  [Winant  el  al.,  1987],  it  is  not  ob¬ 
vious  that  this  balance  will  hold  over  the  shallow  inner  shelf 
where  no  part  of  the  water  column  is  distant  from  the  surface 
and  bottom.  The  moored  array  observations  provide  five  lo¬ 
cations  where  the  two  terms  in  this  balance  can  be  compared 
(Table  6).  In  each  case,  dp/dx  is  estimated  as  the  difference 
between  densities  at  the  same  depth  on  adjacent  moorings 
divided  by  the  mooring  separation.  (Density  at  4.4  m  depth 
for  the  13-m  site  was  computed  by  linearly  interpolating  be- 


LENTZ  ETAL.:  MOMENTUM  BALANCES  ON  THE  NORTH  CAROLINA  INNER  SHELF 


11 


b)  xsy/p0h  and  xby/p0h 


12  17  22  27  1  6  11  16  21  26  1  6  11  16  21  26  31 

Aug  Sep  Oct 


Figure  7.  Time  series  of  largest  terms  in  the  along-shelf  momentum  balance  at  the  13-m  site,  (a)  acceleration  dv/dt  and 
Coriolis  force  fu,  (b)  wind  stress  tsv / p0h  and  bottom  stress  rbv / pah,  and  (c)  pressure  gradient  p^dP/dy  (along  the  5-m 
isobath,  see  Appendix  A.l). 


tween  the  1.5  and  7.6  m  densities.  Figure  2.)  Vertical  shears 
dv/dz ,  estimated  from  current  meters  vertically  bracketing 
the  dp/dx  at  both  moorings,  are  substantial,  corresponding 
to  along-shelf  velocity  differences  of  15  -  40  cm/s  across  the 
water  column. 

The  vertical  shear  in  the  middle  of  the  water  column  is 
approximately  in  thermal  wind  balance  (Table  6,  Figure  13). 
Standard  deviations  of  the  two  terms  in  (6)  are  roughly  equal 
at  each  location  and  increase  toward  the  coast,  consistent 
with  strong  cross-shelf  density  gradients  in  the  vicinity  of 
the  8-m  and  13-m  sites  owing  to  the  Chesapeake  Bay  plume 
and,  during  August,  to  upwelling/downwelling  of  the  pycn- 
ocline.  The  two  terms  in  (6)  are  correlated  at  all  sites.  In¬ 
accuracies  in  the  interpolation  of  density  at  the  13-m  site 
may  contribute  to  the  lower  correlations  for  the  8-m  and  13- 
m  pair.  The  large  separation  between  the  21-m  and  26-m 
sites  (10  km)  relative  to  the  cross-shelf  scale  of  the  plume 
(estimated  from  shipboard  surveys)  probably  contributes  to 


the  lower  correlations  for  this  pair.  Despite  the  difficulties 
of  comparing  vertical  shears  estimated  at  pairs  of  mooring 
sites  with  finite  difference  estimates  of  the  density  gradient 
between  mooring  sites,  the  observations  indicate  that  much 
of  the  subtidal  vertical  shear  variability  is  in  thermal  wind 
balance  in  depths  as  shallow  as  10  m. 

4.  Discussion 

The  momentum  balances  suggest  three  distinct  dynam¬ 
ical  regions:  the  surf  zone,  inner  shelf,  and  mid  shelf.  In 
the  surf  zone,  forcing  by  wave-radiation  stress  gradients 
dominates  both  along-shelf  and  cross-shelf  momentum  bal¬ 
ances.  Over  the  inner  shelf,  offshore  of  the  surf  zone  to 
about  the  13-m  isobath,  the  along-shelf  momentum  bal¬ 
ance  is  primarily  between  surface  wind  and  bottom  stresses, 
with  along-shelf  pressure  gradients  usually  a  smaller  con¬ 
tribution.  The  cross-shelf  momentum  balance  is  between 
cross-shelf  pressure  gradients,  wave-radiation  stress  gradi- 
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Figure  8.  Time  series  of  the  along-shelf  forcing  (rsy / (p0h)  —  ( dP/dy ) / p0)  and  the  response  (dv/dt  +  Tbv /  ( pah ))  at  each 
mooring  site.  The  bottom  stress  rbv  is  estimated  using  the  r  value  seawards  of  the  surf  zone  (9  and  10),  and  dSxy/dx  is  not 
included  in  the  forcing.  Thus,  when  a  site  is  within  the  surf  zone  (  shaded  regions),  discrepancies  may  be  large. 
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b)  (0Sxi»x)/p„h  and  (fv+tsx/p„h) 


Figure  9.  Time  series  of  largest  terms  in  the  cross-shelf  momentum  balance  at  the  8-m  site,  (a)  Coriolis  force  fv  and  wind 
stress  tsx  j p0h  and  (b)  radiation  stress  gradient  ( p0h)~ldSxx/dx  and  (fv  +  tsx / p0h).  Note  the  vertical  scales  differ  in  (a) 
and  (b). 


ents,  Coriolis  forces,  and  to  a  lesser  extent  cross-shelf  wind 
stresses.  At  the  mid  shelf  21-m  and  26-m  sites,  the  along- 
shelf  balance  is  between  wind  stresses,  along-shelf  pressure 
gradients,  bottom  stresses,  and  accelerations.  The  cross¬ 
shelf  momentum  balance  is  geostrophic.  These  simplified 
momentum  balances  are  used  here  to  examine  the  relation¬ 
ship  between  the  forcing  (winds,  waves,  and  along-shelf 
pressure  gradients)  and  the  response  (along-shelf  currents 
and  cross-shelf  pressure  gradients),  and  to  derive  simple 
scalings  relevant  to  other  locations.  The  along-shelf  pressure 
gradient  dP/dy  is  treated  as  a  forcing  rather  than  as  a  re¬ 
sponse  because  the  observed  dP/ dy  (along  both  the  5-m  and 
21-m  isobaths)  are  not  correlated  with  the  local  wind  stress. 
As  indicated  in  Section  3.1,  part  of  the  variability  in  dP/dy 
along  the  5-m  isobath  is  associated  with  the  Chesapeake 
Bay  plume  events  [Rennie  el  al.,  1999],  and  some  of  the 
variability  is  likely  also  a  response  to  the  large-scale  wind 
field  [Wang,  1979;  Noble  and  Butman,  1979;  Yankovsky  and 
Garvine,  1998],  In  either  case,  dP/dy  is  driven  by  large 


scale  processes  not  included  in  the  present  forcing  terms. 

4.1.  Along-shelf  velocity  response 

The  depth-averaged  along-shelf  momentum  balance  (2)  is 
recast  here  to  relate  the  depth-averaged  along-shelf  velocity 
to  the  forcing 

dv  +  rv  ^  1  dP  +  tsv  1  dSxy 
dt  h  p0  dy  p0h  pah  dx 

It  has  been  assumed  that:  turbulent  Reynold's  stresses,  non¬ 
linear  advective  terms,  and  the  Coriolis  term  fu  (see  Section 
3.1)  are  small;  a  linear  drag  law  ((5),  Appendix  A. 3)  pro¬ 
vides  an  accurate  estimate  of  the  bottom  stress  ((/90/i)_1t63/ 
in  (2));  and  the  near-bottom  along-shelf  velocity  Vt,  is  ap¬ 
proximately  equal  to,  or  at  least  well  correlated  with,  the 
depth-averaged  along-shelf  velocity  0.  At  each  site  Vf,  and 
v  are  well  correlated  (0.86  to  0.99)  and  V/,  ranges  from  0.50 
to  0.80  based  on  a  linear  regression  analysis.  Given  that  V[, 
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b)  (dSx)(/3x)/p,h  and  (fv+Tsx/p,h) 


Aug  Sep  Oct 


Figure  10.  Time  series  of  largest  terms  in  the  cross-shelf  momentum  balance  at  the  13-m  site,  (a)  Coriolis  force  fv  and  wind 
stress  tsx / pah,  (b)  radiation  stress  gradient  ( pah)~ldSxx  /  dx  and  (fv  +  tsx  /  p0h),  and  (c)  pressure  gradient  pa~ldP / dx 
and  (fv  +  tsx / p0h  —  (p0h)~ldSxx/  dx).  Correlations  and  regression  coefficients  are  given  in  Table  5. 


and  v  are  well  correlated  the  difference  in  magnitude  may  be 
accounted  for  by  adjusting  r. 

Integrating  (7)  in  time  yields  [Lentz  and  Winant,  1986] 

~  _  [*(  1  dP^T°y  i  dsxy \ 

Vp  J to  V  Po  9y  pah  pah  dx  ) 6 

+  vae  Tf  (8) 

where  Tf  =  h/r  is  a  frictional  time  scale  and  va  =  v(t  = 
t0).  Based  on  (7)  or  (8),  two  factors  determine  the  character 
of  v  and  its  cross-shelf  structure,  the  relative  magnitudes  of 
the  forcing  terms  and  the  forcing  time  scale  relative  to  the 
frictional  time  scale  Tf.  As  noted  above,  dSxy/dx  domi¬ 
nates  the  along-shelf  forcing  in  the  surf  zone.  Seaward  of 
the  surf  zone,  dSxy/dx  is  small  and  the  relative  importance 
of  rsy /h  and  dP/dy  depends  on  their  cross-shelf  structure. 
In  the  absence  of  cross-shelf  variations,  dP/dy  will  become 
increasingly  more  important  relative  to  rsy /h  as  the  depth 
increases  [e.g.  Lentz  and  Winant,  1986],  and  this  is  qualita¬ 


tively  true  in  the  present  observations  (compare  the  cross¬ 
shelf  variation  of  rsy / h  with  that  of  dPj dy  in  Table  2). 

If  the  forcing  varies  on  time  scales  that  are  long  com¬ 
pared  with  Tf  then  the  response  is  basically  frictional.  The 
along-shelf  flow  is  approximately  in  phase  with  the  forc¬ 
ing  and  quasi-steady  in  the  sense  that  the  accelerations 
are  dynamically  negligible.  In  this  case,  from  (7),  vp  ~ 
Tf(-p-'dP/dy  +  rsy  /  p0h  +  (pah)~1dSxy/dx).  In  con¬ 
trast,  if  the  forcing  time  scale  is  short  compared  to  Tf  then 
the  response  lags  the  forcing  and  is  weaker  than  the  steady 
response.  The  three  forcing  terms  (sampled  hourly)  have 
decorrelation  time  scales  of  1  -  2  days.  The  results  in  Sec¬ 
tion  3  suggest  that  roughly 

r  =  5  x  10“3m/s  h  <  3.0 Hsig  (9) 

(within  the  surf  zone), 

5  x  10_4m/s  h  >  3.0-f/gig 


(10) 
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Mooring 

Sites 

Meters  Below 

Surface 

Standard  Deviations 

—g  vp  Vv 

Po  f  dx  dz 

Correlation 

Coefficient 

Days 

8  -  13 

4.4  m 

3.6 

3.3 

0.65 

49 

13-21 

7.6  m 

2.5 

2.5 

0.89 

31 

13-21 

12.7  m 

1.1 

0.9 

0.72 

53 

21  -26 

7.6  m 

0.9 

1.1 

0.62 

65 

21  -26 

12.7  m 

0.8 

0.6 

0.53 

62 

Table  6.  Comparison  of  terms  in  the  thermal  wind  balance  (6). 
significantly  different  from  zero  at  the  95%  confidence  level. 


fv  and  (3P/3x)/p0  26-m  site 


Figure  11.  Time  series  of  largest  terms  in  the  cross-shelf 
momentum  balance  at  the  26-m  site,  fv  and  p~18P/dx. 
Correlations  and  regression  coefficients  between  fv  and 
p~ 1  dP / 8x  are  given  in  Table  5. 


Figure  12.  Time  series  of  the  barotropic  and  baroclinic  con¬ 
tributions  to  the  depth-averaged  cross-shelf  pressure  gradi¬ 
ent  (see  (3))  at  the  13-m  site. 


Standard  deviation  units  are  10  2s  1.  All  correlations  are 


(seawards  of  the  surf  zone). 

Within  the  surf  zone,  assuming  h  =  8  m  or  less,  Tf  is  30 
minutes  or  less,  short  compared  with  the  forcing  time  scale 
and  consistent  with  previous  observations  in  which  surf  zone 
flows  responded  within  a  few  hours  to  changes  in  forcing 
[Feddersen  et  al.,  1998].  Seawards  of  the  surf  zone,  h  =  4  — 
26  m  and  (10)  implies  Tf  =  2  —  14  hours.  Thus,  at  mid  shelf 
Tf  approaches  the  forcing  decorrelation  time  scale  and  the 
flow  response  may  detectably  lag  the  forcing.  The  observed 
time  lags  for  maximum  correlation  between  subtidal  flow  v 
and  forcing  (— dP/dy+rsy /h+h~1dSxy/dx)  are  2  hr,  7  hr, 
and  9-11  hr  for  the  4-m,  8-m,  and  deeper  sites,  respectively. 
The  increase  in  the  magnitude  of  the  lags  with  increasing 
water  depth  is  consistent  with  (7)  and  (9  and  10). 

Estimates  of  vp  from  the  forcing  were  made  using  (8)  and 
the  prescription  for  r  in  (9  and  10).  The  frictional  time  scale 
Tf  in  the  surf  zone  is  shorter  than  the  hourly  sample  rate, 
therefore  Vp  was  set  equal  to  the  forcing  times  Tf  when  the 
4-m  or  8-m  sites  were  in  the  surf  zone.  Unfiltered  hourly 
time  series  of  Vp  and  observed  v  for  each  site  are  similar 
(Figure  14).  Root-mean-square  (rms)  differences  are  about 
10  cm/sec,  correlations  are  0.74  to  0.86,  and  linear  regres¬ 
sion  slopes  range  from  0.92  to  1.20. 

The  vp  estimates  reproduce  the  observed  cross-shelf  vari¬ 
ations  in  v.  For  example,  on  August  20  and  23,  strong  along- 
shelf  currents  are  driven  in  part  by  an  along-shelf  pressure 
gradient  associated  with  the  Chesapeake  Bay  plume  (Sec¬ 
tion  3.1).  The  observed  and  predicted  currents  are  maxi¬ 
mum  at  the  13-m  site  (about  -40  cm/s  and  -70  cm/s  on  20 
and  23  August,  respectively).  Flows  are  weaker  at  the  21- 
m  and  26-m  sites  because  the  pressure  gradient  associated 
with  the  plume  does  not  extend  offshore  to  these  sites.  The 
pressure  gradient  is  present  onshore  of  the  13-m  site,  but  the 
response  decreases  as  the  depth  decreases  from  13-m  to  4- 
m  because  the  pressure  gradient  body  force  is  balanced  by 
bottom  stress  (e.g.  rvp  ~  h  dP/dy  from  (7)).  The  es¬ 
timated  velocity  vp  does  not  reproduce  v  as  well  during  a 
similar  plume  event  around  September  20,  probably  because 
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Figure  13.  Time  series  of  the  terms  in  the  thermal  wind  balance  (6)  (a)  between  the  8-m  and  13-m  sites  (4.4  m  below  surface) 
and  (b)  between  the  13-m  and  21-m  sites  (12.7  m  below  surface).  Correlations  are  given  in  Table  6. 


the  array  does  not  resolve  cross-shelf  variations  in  dP/dy 
during  this  event.  The  vp  estimates  using  the  bottom  stress 
formulation  (9  and  10)  does  reproduce  the  flow  reversal  on 
September  22  (when  dSxy/dx  opposes  the  wind)  between 
the  surf  zone  (v  =  +80  cm/s  at  the  4-m  site)  and  the  inner 
shelf  ( v  =  —50  cm/s  at  the  13-m  site).  The  flow  reversal 
observed  between  the  4-m  and  21-m  sites  on  October  14- 
15  is  also  predicted.  The  crude  drag  formulation  (9  and  10) 
contributes  to  errors  in  the  magnitude  of  vp  in  the  transition 
region  between  the  surf  zone  and  the  inner  shelf.  Errors  in 
the  sign  of  vp  (e.g.  at  the  8-m  site  on  22  September  where 
vp  =  +15  cm/s  and  v  =  —30  cm/s)  may  occur  because 
of  inaccurate  estimation  of  the  small  residual  stress  when 
opposing  stresses  from  wave  breaking  and  wind  are  about 
equal. 

Overall,  these  comparisons  (e.g.  Figure  14)  suggest  that 
(8)  estimates  v  well,  given  the  forcing  and  the  prescription  (9 
and  10)  for  the  bottom  drag  coefficient.  It  is  somewhat  sur¬ 
prising  that  the  linear  drag  formulation  (9  and  10)  works  this 
well,  given  its  likely  dependence  on  waves,  bottom  rough¬ 
ness,  and  stratification.  To  determine  the  sensitivity  of  vp  to 
variations  in  r  (seawards  of  the  surf  zone),  the  rms  difference 
between  vp  and  v  was  computed  as  a  function  of  r  over  the 
range  1  —  20  x  10-4  m/s.  Seaward  of  the  surf  zone,  the  rms 
error  is  not  sensitive  to  r  between  3  —  6  x  10_4m/s.  The  rms 
difference  is  more  sensitive  to  low  values  of  r  than  to  high 


values.  Clearly,  bottom  stress  remains  a  poorly  understood 
aspect  of  the  surf  zone  and  inner  shelf  dynamics. 


4.2.  Cross-shelf  pressure  gradient  response 

The  cross-shelf  pressure  gradient  balances  wave,  wind, 
and  Coriolis  (associated  with  the  along-shelf  flow)  forces 


dP 

dx 


1  dSxx 
h  dx 


+ 


+  Pofv. 


(ID 


Over  mid  shelf  the  present  (Section  3.2)  and  previous  re¬ 
sults  indicate  that  the  cross-shelf  balance  (11)  is  approxi¬ 
mately  geostrophic.  However,  as  the  depth  decreases  the 
cross-shelf  wind  and  radiation  stresses  become  more  impor¬ 
tant.  The  cross-shelf  gradient  of  the  depth-averaged  pressure 
may  be  estimated  from  (11)  given  cross-shelf  distributions 
of  depth,  wave  radiation  stress  Sxx ,  cross-shelf  wind  stress, 
and  the  along-shelf  current  v  (which  may  be  estimated  us¬ 
ing  (7)  and  (8),  given  the  wave  radiation  stress  Sxy ,  along- 
shelf  wind  stress,  and  along-shelf  pressure  gradient).  If,  on 
the  inner  shelf  (seaward  of  the  surf  zone,  where  dSxy/dx 
is  negligible),  the  wind  stress  is  approximately  equal  to  the 
bottom  stress  in  the  along-shelf  momentum  balance  (e.g. 
v  —  rsy / ( p0  r )  in  (7),  see  Section  3.1)  then  the  cross-shelf 
pressure  gradient  (11)  depends  only  on  the  wave  radiation 


cm/s  cm/s  cm/s  cm/s 
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50  8-m  site 


Figure  14.  Time  series  (unfiltered  hourly)  of  depth-averaged  along-shelf  velocities  at  each  site.  Red  curves  are  observed  v 
and  blue  curves  are  vp  estimated  from  forcing  using  (8)  and  (9  and  10).  The  vertical  scale  is  different  for  the  4-m  site. 
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stress  Sxx  and  the  local  wind  stress 

dP  1  dSxx  tsx  frsy 

"o  —  —  r- u - * — 1 — I - •  (12) 

ox  h  ox  h  r 

To  examine  the  relative  magnitude  of  the  three  terms  on 
the  right-hand  side  of  12),  consider  first  the  ratio  of  the  cross¬ 
shelf  wind  stress  to  the  Coriolis  force  is 

rsx/h  _  r  tan (</>) 
frsy /r  fh 

where  an  along-shelf  wind  stress  direction  corresponds  to 
cf>  =  0.  The  depth  h  where  the  magnitude  of  the  two  terms 
are  equal  as  a  function  of  the  wind  sttess  orientation,  for  mid 
latitudes  (  /  =  10-4  s_1)  and  a  typical  drag  coefficient  sea¬ 
ward  of  the  surf  zone  (r  =  5  x  10_4m/s)  are  shown  in 
Figure  15a.  For  winds  oriented  45°  relative  to  the  coastline, 
as  is  often  the  case  at  Duck,  the  cross-shelf  wind  stress  is 
about  half  the  Coriolis  force  in  12-m  depth  and  the  terms  are 
approximately  equal  in  6-m  depth  (consistent  with  observa¬ 
tions,  Figures  10a  and  9a).  The  cross-shelf  wind  stress  will 
be  relatively  more  important  at  lower  latitude  or  if  the  drag 
coefficient  is  larger. 

Observational  studies  concluding  that  the  cross-shelf  mo¬ 
mentum  balance  is  primarily  geostrophic  are  often  in  depths 
greater  than  20  m,  and  do  not  consider  the  contribution  of 
wave -driven  set-down  associated  with  nonbreaking  surface 
waves  in  variable  depth.  In  the  present  observations  the 
wave -radiation  stress  gradients  are  at  least  as  large  as  the 
cross-shelf  wind  stress  and  the  Coriolis  force  in  0(10  m) 
depth  seawards  of  the  surf  zone  (Figures  9b  and  10b).  The 
generality  of  this  result  is  assessed  by  estimating  the  size  of 
h~x  dSxx/dx  relative  to  the  Coriolis  term  frsy  /r  and  the 
cross-shelf  wind  stress  tsx  /h. 

Assuming  normally  incident  waves  for  simplicity,  the  ra¬ 
diation  stress  given  by  (A4)  is 

Sxx  =  E^-^j  (14) 


where  the  phase  speed  c  =  w/ k  and  the  group  velocity  cg  = 
dto/dk  follow  from  the  linear  wave  dispersion  equation 


io2h 

9 


kh  tanh(fc/i) 


where  u>  and  k  are  the  wave  radian  frequency  and  wavenum¬ 
ber,  respectively.  Energy  conservation  (for  the  present  case 
of  nonbreaking  waves)  yields 


OO  Po9-^o o,sig 

Cg  16(tanh(fc/i)  +  khsech2(kh)) 


Figure  15.  (a)  Depth  where  Coriolis  force  associated 

with  wind-driven  along-shelf  flow  \fv\  equals  the  cross¬ 
shelf  wind  sUess  \rsx\/p0h,  as  a  function  of  wind  direc¬ 
tion  (13).  An  along-shelf  directed  wind  corresponds  to 
0°.  (b)  Depth  where  \fv\  equals  the  radiation  stress  gra¬ 
dient  (p0h)~1\dSxx/dx\,  as  a  function  of  H2is  oohx/rsx 
(16).  (c)  Depth  where  |rsxj  equals  \dSxx/dx\,  as  a  func¬ 
tion  of  77/ g  00hx/rsx  (17).  Each  curve  in  b)  and  c)  corre¬ 
sponds  to  the  indicated  wave  period.  Terms  estimated  with 
/  =  10_4s_1,  r  =  5  x  10_4m/s,  and  p0  =  1023  kg/m3. 
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with  HoO’Sig  the  deep  water  significant  wave  height.  Substi¬ 
tution  into  (14)  yields 


^xx 


PoQH^jJjs) 

16 


(15) 


where  s  =  kh  and 


/(*) 


1/2  +  s  sech(s)csch(s) 
tanh(s)  +  s  sech2(s) 


Using  the  dispersion  equation,  s  depends  on  the  nondimen- 
sional  depth  uj2h/  g.  Differentiation  yields 


1  dSxx 
h  dx 


_ n  tt2 

16^°  °°,sis 


(16) 


where 


«(*) 


/'(«) 


coth(s) 

s(tanh(s)  +  ssech2(s)) 


The  ratio  of  the  radiation  stress  term  to  the  Coriolis  term  is 
thus 


iasxx\  (fT°y 

h  dx  )  \  r 


Po  T 

16  gf 


tj2 

-°sig,oo 


h  x 


rsv 


■u)4  q(s).  (17) 


Similar  calculations  for  the  ratio  of  dSxx/dx  and  rs 
yield 


dSxx  sx  =  Po  _  H,^hx 

•  1  (2.  si-  agr.  *  '  ' 


dx 


16 


(18) 


where 


p(s)  = 


m 

tanh(s)  +  s  sech2(s) 


The  functions  q  and  p  depend  relatively  strongly  on  uj2h/g 
(Figure  16).  For  example,  in  shallow  water  p  ~  ft,-5/2  and 
q  ~  ft-3/2.  In  contrast,  for  wind  stresses  that  do  not  vary 
in  the  cross-shelf  direction,  the  cross-shelf  wind  stress  term 
tsx  / h  varies  as  h~ 1  and  the  Coriolis  term  frsy  /r  is  con¬ 
stant.  These  ratios  ((17)  and  (18))  suggest  wave  radiation 
stresses  will  dominate  in  very  shallow  water,  even  with  non¬ 
breaking  waves. 

The  depth  where  the  ratio  of  the  radiation  stress  to  the 
Coriolis  term  equals  unity  (from  (17))  is  shown  in  Fig¬ 
ure  15b,  for  mid  latitudes  and  a  typical  value  of  the  drag  co¬ 
efficient  (r  =  5  x  10-4  m/s),  as  a  function  of  H2ig  oohx / rsy 
for  several  wave  periods.  At  Duck,  fT2ig  ^  and  rsy  are  corre¬ 
lated  such  that  their  ratio  is  typically  0(20  —  40)  m4/N.  The 
beach  slope  between  the  4-m  and  13-m  sites  is  roughly  0.01, 
so  H2lg  (yohx / rsy  ~  0.2  -  0.4.  Wave  periods  are  0(10  s), 
so  the  ratio  (17)  is  unity  in  about  15-20  m  depth  (Figure 


Figure  16.  The  functions  q(  16)  and  p  (18)  versus  the  nondi- 
mensional  depth  u)2h/g. 


15b),  qualitatively  consistent  with  the  observation  that  the 
terms  are  about  equal  in  13-m  depth  (Figure  10b).  For 
fixed  i72ig  00hx/Tsv ,  the  radiation  stress  term  is  important 
in  deeper  water  as  the  surface  wave  period  (and  hence  wave¬ 
length)  increases  (Figure  15b),  reflecting  the  dependence  of 
q  on  the  nondimensional  depth  (Figure  16). 

At  mid  latitudes,  when  the  wind  direction  is  not  close  to 
cross-shelf  directed  (e.g.  when  the  wind  direction  <  60°  in 
Figure  15a),  the  cross-shelf  wind  stress  exceeds  the  Coriolis 
term  only  when  h  <  10  m.  However,  the  Sxx  gradient  term 
also  usually  exceeds  the  Coriolis  term  in  depths  less  than  10 
m  (Figure  15b)  and  the  Sxx  gradient  term  increases  rapidly 
as  the  depth  decreases  further  (Figure  16).  Under  these  cir¬ 
cumstances,  the  Coriolis  force  dominates  on  the  mid-shelf, 
the  Coriolis  and  radiation  stress  terms  are  both  important 
on  the  inner  shelf,  and  the  cross-shelf  wind  stress  is  never 
a  dominant  term.  In  the  contrasting  situation  of  a  cross¬ 
shelf  wind,  the  Coriolis  term  vanishes  and  the  cross-shelf 
wind  stress  is  larger  than  the  radiation  stress  term  offshore 
(Figure  15c).  The  depth  where  tsx /h  equals  the  radiation 
stress  term  depends  on  77 2 g  vJix  /tsx  and  the  wave  period 
((18)  and  Figure  15c).  The  wave  radiation  stress  gradient 
from  an  energetic  ocean  swell  (e.g.  15  s  period  waves  with 
-Hsig,oo  =  3m)  shoaling  on  a  moderately  sloping  shelf  (e.g. 
hx  =  0.005)  will  exceed  the  force  of  a  strong  cross-shelf 
wind  stress  (0.4  N/m2,  10  m  wind  speed  15  m/s)  in  less 
than  30  m  water  depth  (Figure  15c). 

The  above  results  (e.g.  Figure  15b  and  15c)  suggest  that 
dSxx/dx  ,  tsx /h,  and  frsy /r  may  be  similar  in  magni¬ 
tude  seawards  of  the  surf  zone,  in  depths  of  O(10-30m). 
Within  the  surf  zone,  tsx /h  and  frsy  /r  are  negligible  be¬ 
cause  dSxx  / dx  is  orders  of  magnitude  larger  than  seaward 
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Figure  17.  Summary  schematic  of  the  dominant  cross-shelf 
and  along-shelf  momentum  balances  in  the  surf  zone,  inner 
shelf,  and  mid  shelf  regions.  Dashed  lines  indicate  transition 
zones  where  momentum  balances  may  be  more  complicated. 

of  the  surf  zone  [Longuet-Higgins  and  Stewart,  1964;  Bowen 
el  al.,  1968],  Although,  dSxy/dx  was  neglected  offshore  of 
the  surf  zone  12),  it  can  drive  substantial  along-shelf  currents 
in  the  surf  zone.  Hence,  in  the  surf  zone  (11)  reduces  to 

9P _ 1  8SXX  f  8Sxy 

dx  h  dx  r  dx 

using  (7)  with  a  surf  zone  balance  between  radiation  stress 
gradients  and  bottom  friction.  The  ratio  of  the  two  terms  on 
the  right-hand  side  of  (19)  is 

(. f/r)dSxy/dx  =  fh  sin(2  9) 

( dSxx/dx)/h  r  (cos(2$)  +  2) 

using  (A2)  and  (A4),  assuming  shallow  water  waves  (c  ~ 
cg)  and  similar  cross-shelf  scales  for  the  divergences,  and 
with  6  =  0  corresponding  to  normally  incident  waves. 

The  ratio  of  the  frictional  and  Coriolis  time  scales  (fh/r  = 
fTf )  is  about  0.1  for  h  =  5  m,  r  =  5  x  10-3  m/s,  and 
/  =  10-4  1/s.  The  other  fraction  involving  the  incident 
wave  angle  9  is  always  less  than  0.6,  and  for  waves  within 
10°  of  normally  incident  is  0.1  or  less.  Thus,  (19)  suggests 
that  under  most  surf  zone  conditions  the  pressure  gradient 
associated  with  wave -driven  set-up  will  be  at  least  an  order 
of  magnitude  larger  than  the  geostrophic  pressure  gradient 
associated  with  the  wave-driven  along-shelf  current  jet. 

5.  Summary 

Estimates  of  terms  in  the  along-shelf  and  cross-shelf 
momentum  balances  indicate  three  dynamically  distinct  re¬ 
gions:  the  surf  zone,  the  inner  shelf  between  the  offshore 
edge  of  the  surf  zone  and  the  13-m  isobath,  and  the  mid 
shelf  extending  offshore  of  the  13-m  isobath  (Figure  17). 
Consistent  with  previous  studies,  breaking  surface  gravity 
waves  provide  the  dominant  forcing  in  the  surf  zone.  The 
cross-shelf  divergence  in  the  cross-shelf  component  of  the 


wave -radiation  stress  dSxx  /  dx  is  much  larger  than  the  other 
estimated  terms  suggesting  it  is  balanced  by  a  cross-shelf 
pressure  gradient  (i.e.  wave  setup  that  could  not  be  esti¬ 
mated  with  these  data).  The  present  estimates  also  support 
previous  conclusions  [ Thornton  and  Guza,  1986;  Feddersen 
el  al.,  1998]  that  the  cross-shelf  divergence  in  the  along-shelf 
component  of  the  wave-radiation  stress  dSxy/dx  is  largely 
balanced  by  bottom  stress.  This  balance  requires  a  linear 
drag  coefficient  for  the  subtidal  flow  in  the  surf  zone  that  is 
about  10  times  larger  than  the  drag  coefficient  seaward  of  the 
surf  zone. 

The  cross-shelf  momentum  balance  at  mid  shelf  is  pre¬ 
dominantly  geostrophic,  the  Coriolis  force  due  to  the  along- 
shelf  current  balances  the  cross-shelf  pressure  gradient,  as 
found  elsewhere  [Brown  el  al.,  1985, 1987;  Lee  et  al,  1984, 1989]. 
At  this  site  the  along-shelf  flow  is  geostrophic  into  fairly 
shallow  water,  the  21-m  isobath.  The  along-shelf  momen¬ 
tum  balance  at  mid  shelf  is  more  complex.  The  wind 
stress  and  the  along-shelf  pressure  gradient  terms  are  sim¬ 
ilar  in  magnitude  and  are  balanced  by  both  accelerations  and 
bottom  stress,  consistent  with  previous  studies  [Allen  and 
Smith,  1981;  Lentz,  and  Winant,  1986;  Lee  etal.,  1984,  1989; 
Lentz,  1994]. 

Between  the  surf  zone  and  mid  shelf  is  a  transition  re¬ 
gion  referred  to  as  the  inner  shelf,  here  extending  roughly 
from  the  4-m  to  the  13-m  isobath.  The  along-shelf  momen¬ 
tum  balance  over  the  inner  shelf  is  predominantly  between 
wind  and  bottom  stresses,  consistent  with  the  few  previous 
studies  of  the  momentum  balance  in  10  m  to  15  m  of  water 
[Lentz  and  Winant,  1986;  Lee  etal.,  1989].  Along-shelf  pres¬ 
sure  gradients  were  important  over  the  inner  shelf  during  a 
few  events  associated  with  a  low-salinity  plume  from  Chesa¬ 
peake  Bay,  approximately  100  km  north  of  the  study  site 
[Rennie  et  al.,  1999].  In  the  cross-shelf  direction,  the  Cori¬ 
olis  force  due  to  the  along-shelf  current,  8Sxx/dx  from  the 
shoaling  of  unbroken  surface  gravity  waves,  and  the  cross¬ 
shelf  pressure  gradient  all  have  similar  magnitudes.  The 
cross-shelf  wind  stress  is  never  a  dominant  term.  The  cross¬ 
shelf  wind  stress  is  similar  in  magnitude  to  the  Coriolis  force 
at  the  8-m  site,  consistent  with  the  results  of  Lee  et  al.  [1989] 
at  10  m  water  depth  in  the  South  Atlantic  Bight.  However 
the  radiation  stress  gradient  dSxx/dx,  a  term  not  considered 
by  Lee  et  al.  [1989],  dominates  at  the  8-m  site.  Assuming 
the  cross-shelf  pressure  gradient  balances  the  other  terms, 
the  pressure  gradient  may  be  thought  of  as  a  superposition 
of  a  geostrophic  along-shelf  flow  and  wave  setdown.  The 
strongest  along-shelf  currents  coincided  with  the  largest  sur¬ 
face  waves,  and  geostrophic  setup  (southeastward  flows)  and 
wave  setdown  approximately  cancelled  in  13-m  depth,  re¬ 
sulting  in  a  nearly  flat  mean  sea  surface  during  some  events. 

On  both  the  inner  and  mid  shelf,  the  depth-averaged 
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cross-shelf  pressure  gradient  includes  a  baroclinic  compo¬ 
nent  (cross-shelf  density  gradient)  that  balances  about  one- 
third  to  one-half  of  the  barotropic  component  (cross-shelf 
sea-surface  slope).  This  contrasts  with  results  for  winter  in 
the  South  Atlantic  Bight  where  of  estimates  of  the  baroclinic 
component  (from  shipboard  CTD  surveys)  were  small  com¬ 
pared  to  the  barotropic  component  [Lee  et  al.,  1984,  1989]. 
The  baroclinic  cross-shelf  pressure  gradient  is  in  thermal 
wind  balance  with  the  vertical  shear  in  the  along-shelf  flow 
in  water  as  shallow  as  10  m. 

The  depth-averaged  cross-shelf  velocities  are  too  small 
to  estimate  accurately  from  the  observations,  but  the  ten¬ 
dency  for  the  along-shelf  momentum  balance  to  close  with¬ 
out  including  the  Coriolis  force  fu  (Figures  8  and  14)  sug¬ 
gests  fu  is  not  a  large  term  in  the  along-shelf  momen¬ 
tum  balance.  Several  previous  studies  have  found  simi¬ 
lar  results  for  the  mid  and  inner  shelf  [Allen  and  Smith , 
1981;  Lentz  and  Winant,  1986;  Lee  et  al.,  1989],  though 
Lee  et  al.  [1984,  1989]  find  fu  is  large  over  the  outer  shelf 
in  the  South  Atlantic  Bight.  In  general,  the  depth-averaged 
cross-shelf  flow  remains  poorly  understood.  If  fu  is  small, 
the  along-shelf  and  cross-shelf  depth-averaged  momentum 
balances  decouple  and  the  along-shelf  momentum  balance 
in  these  three  regions  provides  a  simple  estimate  of  the 
depth-averaged  along-shelf  current  in  terms  of  the  forcing, 
i.e.  along-shelf  wind  and  wave  stresses  and  the  along-shelf 
pressure  gradient.  This  simple  estimate  reproduces  accu¬ 
rately  the  observed  depth-averaged  along-shelf  current  and 
its  cross-shelf  variation  from  the  surf  zone  to  mid  shelf  given 
a  linear  drag  coefficient  of  5  x  10-3  m/s  within  the  surf  zone 
and  5  x  10-4  m/s  seaward  of  the  surf  zone. 

Appendix 

The  procedure  for  estimating  terms  in  the  depth-averaged 
momentum  balances  (1)  and  (2)  is  described  in  Appendix 
A.l.  Uncertainties  in  the  estimates  are  discussed  in  Ap¬ 
pendix  A. 2  and  bottom  stress  estimates  from  the  linear  drag 
law  are  compared  to  log-profile  estimates  in  Appendix  A. 3. 

Al.  Estimation  of  Terms 

Terms  in  the  momentum  balances  are  estimated  using  the 
hourly  data  and  then  low-pass  filtered  (half-power  point  38 
hours)  to  focus  on  subtidal  variability.  Vertical  integrals  are 
estimated  using  a  trapezoidal  rule  and  assuming  no  vertical 
variations  near  the  boundaries  to  extrapolate  to  the  surface 
and  bottom. 

Accelerations 

Time  derivatives  are  estimated  as  centered  differences 
over  two  hour  intervals. 

Pressure  Gradients 


The  depth-average  pressure  gradient  in  (1,2)  is  estimated 
from  the  bottom  pressure  and  density  data  by  substituting 
(4)  into  (3)  or  the  equivalent  equation  for  dP/dy.  The  den¬ 
sity  term  in  (3)  is  estimated  more  easily  from  the  moored 
observations  by  noting  that 


dp  z.  , 

9s(1  +  h)dz 


(Al) 

The  integrals  on  the  right-hand  side  are  estimated  for  each 
mooring  site  where  density  estimates  are  available.  The 
cross-shelf  gradients  of  the  bottom  pressure  and  the  density 
integrals  are  estimated  as  finite  differences  roughly  centered 
on  the  site  of  interest.  For  example,  the  cross-shelf  pressure 
gradient  estimate  for  the  21-m  site  is  the  difference  between 
bottom  pressure  and  density  estimates  from  the  26-m  and 
5-m  sites  (bottom  pressure  at  the  5-m  site  and  density  from 
the  FRF  pier  site)  divided  by  the  separation.  The  13-m  esti¬ 
mate  is  from  the  21-m  and  5-m  sites  and  the  26-m  estimate  is 
from  the  33-m  and  5-m  sites  (Figure  2).  Accurate  cross-shelf 
pressure  gradients  can  be  made  only  for  the  13-m,  21-m  and 
26-m  sites.  The  bottom  pressure  and  density  measurements 
at  the  13-m  and  5-m  sites  were  too  close  together  (separa¬ 
tion  1  km)  to  make  reliable  estimates  of  the  depth-average 
pressure  gradient.  There  are  bottom  pressure,  but  no  density 
measurements,  offshore  of  the  26-m  site  (Figure  2).  There¬ 
fore  the  density  gradient  term  in  (3)  for  the  26-m  site  is  esti¬ 
mated  using  observations  from  the  21-m  and  26-m  sites. 


Along-shelf  pressure  gradients  are  estimated  along  the  5- 
m  and  the  21-m  isobaths  using  data  from  the  along-shelf 
pressure/density  array  (Figure  1).  Density  along  the  5-m 
isobath  is  assumed  vertically  uniform  because  temperature 
and  conductivity  were  measured  at  only  one  depth.  Pressure 
gradients  along  the  5-m  isobath  are  estimated  as  differences 
between  measurements  from  the  sites  17  km  north  and  south 
of  the  central  line.  Pressure  gradients  along  the  21-m  isobath 
are  estimated  as  differences  between  measurements  from  the 
sites  32  km  north  and  27  km  south  of  the  central  line.  Pres¬ 
sure  gradient  estimates  from  different  pairs  of  bottom  pres¬ 
sure  sensors  along  either  the  5-m  or  21-m  isobaths  yielded 
similar  estimates.  The  pressure  gradient  estimate  along  the 
5-m  isobath  is  compared  with  the  other  terms  in  the  along- 
shelf  momentum  balance  from  the  4-m,  8-m  and  13-m  sites 
and  the  pressure  gradient  along  the  21-m  isobath  to  terms 
from  the  21-m  and  26-m  sites. 


Wind,  Bottom,  and  Wave-Radiation  Stresses 

The  wind  stress  is  estimated  using  a  neutral  drag  law 
[Large  and  Pond ,  1981]  and  the  FRF  pier  winds.  Other  neu¬ 
tral  and  non-neutral  bulk  estimates  [Fairall  et  al.,  1996]  give 
nearly  identical  wind  stresses  [Austin,  1998].  No  attempt 
has  been  made  to  account  for  the  influence  of  waves  and  the 
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shallowness  of  the  water  in  the  drag  law  used  because  the 
dependence  of  the  drag  coefficients  on  these  effects  is  not 
understood  well  [ Geernaert ,  1988].  The  FRF  pier  winds  are 
used  because  there  is  a  continuous  time  series  over  the  dura¬ 
tion  of  the  study  and  the  pier  wind  measurements  are  nearly 
identical  to  the  wind  measurements  at  the  21-m  site  (magni¬ 
tude  of  the  vector  correlation  is  0.98). 

Bottom  stress  is  estimated  assuming  a  linear  drag  law  (5) 
with  r  =  5  x  10-4  m/s.  This  choice  is  motivated  by  sim¬ 
plicity  and  the  poor  understanding  of  the  factors  influencing 
bottom  stress  in  the  surf  zone  and  inner  shelf.  Bottom  stress 
estimates  from  (5)  are  compared  with  log-profile  estimates 
using  bottom  tripod  measurements  at  the  21-m  site  in  Ap¬ 
pendix  A. 3. 

Wave -radiation  stresses  are  estimated  in  8-m  depth  from 
the  FRF  wave-directional  array  data  and  a  directional-moment- 
estimation  technique  [Elgar  et  al.,  1994].  To  estimate  Sxy  at 
other  locations,  the  depth  and  surface  wave  properties  are  as¬ 
sumed  not  to  vary  in  the  along-shelf  direction,  and  the  wave 
field  is  assumed  narrow-banded  in  both  frequency  and  direc¬ 
tion.  In  this  case  [Longuet-Higgins  and  Stewart,  1964] 


Sxy  =  Ecg  sin(20)/(2c),  (A2) 

where  E  =  p0gH/ig/ 16  is  the  wave  energy,  Hs ;g  is  the  sig¬ 
nificant  wave  height,  pQ  is  a  reference  density,  g  is  gravita¬ 
tional  acceleration,  cg  and  c  are  the  (linear  theory)  group  and 
phase  velocities,  respectively,  at  the  peak  wave  frequency. 
The  angle  6  in  8-m  depth  is  chosen  such  that  Sxy  from  (A2) 
equals  the  measured  Sxy  in  8-m  depth  (using  values  of  E, 
and  of  c  and  cg  based  on  the  peak  frequency,  measured  in 
8-m  depth)  [ Thornton  and  Guza,  1986]. 

The  cross-shelf  gradient  in  Sxy  at  each  measurement  site, 
using  (A2)  and  Snell’s  Law  is  given  by 

gXV  =  e  sin(fl)/c  (A3) 

where  the  wave  dissipation  e  =  d(Ecg  cos (9))/dx  is  esti¬ 
mated  using  the  model  of  Whitford  and  Thornton  [ Church 
and  Thornton,  1993],  with  D  =  0.72  and  7  =  0.3  [Chen 
et  al.,  1997].  The  modeled  dissipation  rate  depends  linearly 
on  the  local  bottom  slope  (estimated  from  bathymetric  sur¬ 
veys,  Figure  2)  and  on  the  local  wave  energy  (obtained  by 
applying  linear  theory  to  the  bottom  pressure  spectrum  at 
sea-swell  frequencies  measured  near  each  mooring/tower.) 
Cross-shelf  gradients  of 


Sxx  =  E 


^(l+c„s2(#))-l 


(A4) 


were  estimated  similarly  using  estimates  of  the  dissipation 
from  the  model  and  estimates  of  E  and  the  peak  wave  fre¬ 
quency  from  observations  at  each  site.  When  there  is  no 


wave  breaking,  Sxy  gradients  vanish,  but  Sxx  gradients  are 
nonzero  [e.g.  Longuet-Higgins  and  Stewart,  1964,  and  (16)]. 

A2.  Uncertainties 

Uncertainties  in  the  estimates  described  above  cannot  be 
quantified  accurately  because  instrument  performance  in  the 
field  is  not  understood  well  and  because  spatial  scales  of 
variation  are  unknown.  Nevertheless,  a  qualitative  discus¬ 
sion  of  uncertainties  provides  some  context  for  interpreting 
the  estimates  presented  in  Section  3. 

Errors  in  the  depth-averaged  velocity  estimates  include 
current  speed  and  direction  measurement  errors  and  errors  in 
estimating  the  depth-averaged  currents  from  a  few  vertically 
spaced  current  measurements.  The  current  meters  have  a  re¬ 
ported  accuracy  of  2-3  cm/s  [Beardsley,  1987;  Guza  et  al, 
1988],  Uncertainties  in  orientation  have  a  larger  impact  on 
the  cross-shelf  velocities  than  along-shelf  velocities  because 
the  flows  are  strongly  polarized  (Table  1).  The  largest  source 
of  error  in  estimating  depth-averages  is  probably  extrapola¬ 
tion  of  the  current  profiles  to  the  surface  and  bottom.  How¬ 
ever,  comparison  of  the  standard  estimate  at  the  21-m  site 
with  an  estimate  incorporating  near-surface  OSCAR  [Shay 
et  al.,  1998]  and  near  bottom  tripod  (Appendix  A. 3)  current 
measurements  yielded  rms  differences  in  depth-averaged  ve¬ 
locities  of  1  cm/s.  Maximum  differences  were  5  cm/s  for 
hourly  data  and  3  cm/s  for  low-pass  filtered  data.  This  sug¬ 
gests  depth-averaging  at  this  site  does  not  increase  substan¬ 
tially  the  uncertainty  in  the  estimates  over  the  uncertainty 
in  the  individual  current  measurements.  Current  meter  un¬ 
certainties  include  wave -induced  biases  that  may  affect  all 
instruments  on  a  mooring  and  thus  it  can  not  be  assumed 
that  depth-averaging  will  reduce  the  uncertainty  from  that 
for  an  individual  current  meter.  Therefore  Su  =  3  cm/s  is 
the  assumed  uncertainty  in  both  the  currents  and  the  depth- 
averaged  currents.  The  subtidal  accelerations  have  an  es¬ 
timated  uncertainty  of  Su/  At,  where  a  At  of  9  hours  was 
chosen  as  approximately  a  quarter  of  the  38  hour  cut-off  pe¬ 
riod  for  the  low -pass  filter  (Table  Al).  Estimated  uncertainty 
in  the  Coriolis  term  is  fSu. 

Parallel  plate  pressure  ports  were  added  to  the  Seagauges 
to  reduce  flow  noise.  Based  on  flume  tests  (flow  speeds  of  6 
-  40  cm/s),  and  on  ocean  deployments  of  several  sensors  at 
the  same  site,  the  pressure  measurements  have  a  relative  ac¬ 
curacy  of  0.1 -0.2  mb.  The  influence  of  waves  on  the  pressure 
measurements  is  not  known.  Comparisons  of  the  Seagauges 
and  Setra  pressure  sensors  at  the  5-m  and  21-m  central  line 
sites  indicate  relative  accuracies  between  these  pairs  of  ~1 
mb.  Uncertainties  in  the  bottom  pressure  gradient  terms  are 
100(SPb/ Ax) / Po  (100  is  the  conversion  from  mb  to  N /m2), 
where  SPb  is  assumed  to  be  0.2  mb  for  Seagauge  pairs  and  1 
mb  for  Seagauge-Setra  pairs.  Uncertainties  in  the  density  are 
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Term 

4-m 

8-m 

13-m 

21-m 

26-m 

5u/  At 

0.1 

0.1 

0.1 

0.1 

0.1 

fSu 

0.3 

0.3 

0.3 

0.3 

0.3 

5Pb/[p0  Ax) 

- 

11.7 

0.4 

0.7 

0.3 

ghSp/(p0Ax) 

- 

0.9 

0.3 

0.1 

0.1 

6Pb/(p0Ay) 

0.06 

- 

- 

0.03 

- 

ghSp/(p0Ay) 

0.02 

- 

- 

0.08 

- 

5ts  Kp0ti) 

0.7 

0.4 

0.2 

0.1 

0.1 

5rb  !(p  0h) 

0.4 

0.2 

0.1 

0.1 

0.1 

Table  Al.  Al.  Estimated  uncertainties  in  the  terms  of  the 
cross-shelf  and  along-shelf  momentum  balances  (1)  and  (2). 
Units  are  10-5m/s2. 

estimated  to  be  0.1  kg/m3  based  on  comparisons  with  ship¬ 
board  CTD  measurements  and  between  adjacent  instruments 
on  the  same  mooring  when  the  water  column  is  thought  to 
be  well  mixed.  Uncertainties  arise  primarily  from  drift  of  the 
conductivity  cells  [Alessi  et  al.,  1996],  Uncertainties  in  the 
depth-averaged  densities  due  to  having  only  two  instruments 
at  the  northern  and  southern  21-m  sites  were  estimated  by 
comparing  depth-averages  using  all  four  instruments  at  the 
central  21-m  site  with  estimates  using  only  two  instruments. 
RMS  differences  were  0.25  kg/m3.  Thus,  Sp  is  assumed 
to  be  0.1  kg/m3  except  for  the  along-shelf  gradient  at  the 
21-m  site  where  it  is  assumed  to  be  0.25  kg/m3.  Uncer¬ 
tainties  in  the  baroclinic  pressure  gradient  are  estimated  as 
gh8p/(p0Ax).  There  is  additional  uncertainty  associated 
with  cross-shelf  variations  in  dP/dy  between  5  m  and  13 
m  and  between  21m  and  26  m  because  dP/dy  is  only  esti¬ 
mated  along  two  isobaths. 

Hourly  wind  stress  estimates  from  the  FRF  pier  and  the 
21-m  site  (separation  ~5  km)  have  rms  differences  of  2  —  3  x 
10-2  N/m2.  These  differences  are  presumably  from  both 
instrument  inaccuracies  and  spatial  variations  in  the  wind 
which  is  assumed  to  be  spatially  uniform.  Thus,  the  esti¬ 
mated  uncertainty  is  8ts /(p0h)  with  8ts  =  3  x  10-2  N/m2. 
Another  major  source  of  uncertainty  is  the  drag  coefficient 
in  shallow  water  and  its  dependence  on  sea  state  and  wind 
direction  (fetch).  A  summary  of  drag  coefficient  estimates 
in  shallow  water  by  [Geernaert,  1988]  suggests  uncertain¬ 
ties  of  about  20%.  Recent  evidence  suggests  the  drag  coef¬ 
ficient  over  the  inner  shelf  may  be  different  for  onshore  and 
offshore  winds  [Friedrichs  and  Wright,  1998]. 

Uncertainties  in  the  bottom  stress  term  owing  to  uncer¬ 
tainties  in  the  current  measurements  are  rSu/h.  However, 
of  more  concern  is  the  dependence  of  r  on  factors  such  as 
waves,  bottom  roughness,  and  stratification,  and  more  gen¬ 
erally  the  appropriateness  of  a  linear  drag  law  for  parameter¬ 
izing  bottom  stress.  A  comparison  of  the  linear  drag  and  log- 
profile  estimates  of  bottom  stress  at  the  21-m  site  is  given  in 


Figure  Al.  Time  series  of  along-shelf  bottom  stress  at  the 
21-m  site  estimated  using  a  linear  drag  law  with  r  =  5  x 
10_4m/s  and  log-profiles  from  either  a  bottom  tripod  or  a 
bottom  tetrapod.  Tetrapod-based  estimates  are  shown  during 
August,  when  a  tripod  and  tetrapod  were  both  deployed. 

Appendix  A. 3. 

The  radiation  stress  gradient  estimates  are  crude.  Errors 
are  owing  to  uncertainties  in  the  bottom  slope,  errors  in  lin¬ 
ear  theory,  and  errors  in  model  based  estimates  of  the  break¬ 
ing  wave  dissipation  rate.  The  radiation  stress  gradient  es¬ 
timates  for  nonbreaking  waves  are  probably  accurate  within 
50%,  but  may  be  less  accurate  for  breaking  waves.  When 
the  fraction  of  waves  breaking  is  low,  but  not  negligible,  as 
occurs  in  the  region  bordering  the  seaward  edge  of  the  surf 
zone,  the  true  Sxx  gradient  changes  sign  and  the  estimates 
may  have  larger  fractional  errors  and/or  the  wrong  sign. 

A3.  Bottom  Stress  Estimates 

The  linear  drag  formulation  used  to  estimate  bottom  stress 
is  crude  because  it  does  not  account  for  factors  such  as  vari¬ 
ability  in  bottom  roughness  and  surface  gravity  waves.  Bot¬ 
tom  tripod  and  tetrapod  [Kim  et  al.,  1997]  deployments  at 
the  21-m  site  provide  independent  bottom  stress  estimates. 

The  1-m  tall  tetrapods  supported  four  Marsh-McBirney 
electromagnetic  current  meters  sampling  at  1  Hz  for  1024 
seconds  every  4  hours,  and  an  altimeter  that  determined 
the  sensor  elevations  above  the  sea  floor.  There  were  30- 
day  tetrapod  deployments  during  August  and  during  Oc¬ 
tober.  The  5-m  tall  bottom  tripods  supported  five  BASS 
current  meters  [Williams  et  al.,  1987]  and  seven  thermis¬ 
tors.  Currents  were  sampled  nearly  continuously  at  approx¬ 
imately  1.5  Hz.  There  were  30-day  tripod  deployments  dur¬ 
ing  July-August  and  during  September.  There  is  generally 
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Deployment 

Filtered 

Tby  Tbx 

Duration 

Days 

Current  Meter 

Elevation  (m) 

Full  time  series 

0.81* 

0.45* 

74 

Aug  tetrapod 

0.79* 

0.39 

20 

0.06,  0.35,  0.66,  0.95 

Oct  tetrapod 

0.87* 

0.56 

28 

0.05,  0.34,  0.65,  0.94 

Aug  tripod 

0.94* 

0.75 

12 

0.24,  1.20,2.55,  4.44 

Sep  tripod 

0.89* 

0.52 

26 

0.24,  0.60,  1.20,2.55,  4.44 

Table  A2.  Correlations  between  subtidal  bottom  stresses  estimated  at  the  21-m  site  from  a  linear  drag  law  (using  mooring 
observations)  and  log-profiles  (using  tripod  and  tetrapod  observations).  The  duration  of  measurements  used  in  the  correlations 
and  current  meter  elevations  above  the  sea  floor  are  also  given.  *  denotes  correlations  significantly  different  from  zero  at  the 
95%  confidence  level. 


good  agreement  between  the  nearest-bottom  VMCM  (eleva¬ 
tion  1.5  m)  on  the  21-m  mooring  and  the  BASS  tripod  and 
EMCM  tetrapod  current  measurements,  with  the  exception 
of  the  October  tetrapod  deployment  when  along-shelf  tetra¬ 
pod  velocities  exceeded  the  VMCM  velocities,  despite  being 
closer  to  the  bottom.  It  is  unclear  which  measurements  are 
correct. 

Bottom  stresses  are  estimated  from  the  tripod  and  tetra¬ 
pod  data  using  a  log-profile  technique  [e.g.  Kim  el  al.,  1997]. 
Only  the  two  sensors  closest  to  the  sea  floor  are  used  to  esti¬ 
mate  bottom  stress  as  there  is  curvature  in  the  speed  profiles, 
perhaps  because  the  water  column  is  often  stratified  within  a 
few  meters  of  the  bottom.  Using  the  lowest  three,  rather  than 
the  lowest  two  sensors  has  little  effect  on  the  bottom  stress 
estimates  except  during  October  when  the  profile  curvature 
is  persistent  and  strong. 

The  along-shelf  bottom  stress  estimates  from  log-profiles 
and  the  linear  drag  law  agree  well  (Figure  Al).  Zero-lag 
correlations  are  0.79  or  greater  for  each  tripod/tetrapod  de¬ 
ployment  and  for  the  combined  time  series  (Table  A2).  The 
most  notable  discrepancies  are  the  three  events  in  October 
when  the  log-profile  estimates  yield  bottom  stresses  that  are 
about  twice  as  large  as  the  linear  drag  law  (Figure  Al).  It  is 
unclear  which  estimate  is  more  accurate.  If  log-profile  esti¬ 
mates  of  bottom  stress  are  used  instead  of  linear  drag  law  es¬ 
timates,  the  correlation  between  forcing  and  flow  response  at 
the  21-m  site  is  slightly  (though  not  significantly)  increased, 
from  0.88  (Table  3)  to  0.91.  Correlations  between  linear  and 
log-profile  bottom  stress  estimates  are  smaller  for  the  weaker 
cross-shelf  component  of  bottom  stress  (Table  A2),  but  the 
magnitudes  of  both  estimates  are  similar  and  much  smaller 
than  other  terms  in  the  cross-shelf  momentum  balance  (Ta¬ 
ble  4).  As  log-profile  estimates  are  available  only  at  the  21-m 
site,  the  linear  drag  law  is  used  at  all  sites. 
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